Formação e integração empresarial de quadros (FIEQ) 


Por despacho conjunto dos secre- 
tários de Estado da População e 
Emprego e das Indústrias Extrac- 
tivas e Transformadoras foi criado 
o programa experimental FIEQ, a 
ser lançado pelo IAPMEI, com o 
objectivo de desenvolver um esque- 
ma de apoio às PMEs nos domi- 
nios tecnológico, de organização e 
gestão. 

Para melhor elucidação dos pos- 
síveis interessados transcreve-se de 
seguida o citado despacho, de 17 
de Novembro de 1978. 


Despacho conjunto 


1 — Graves carências das peque- 
nas e médias empresas industriais 
(PME's) nos dominios tecnológicos 
de organização e gestão coexistem 
com numerosas situações de de- 
semprego e subemprego de quadros 
técnicos, nomeadamente jovens re- 
cém-licenciados, que, embora sem 
experiência profissional, poderão 
prestar serviços úteis àquelas em- 
presas se lhes for proporcionado 
um estágio que lhes assegure con- 
dições de melhor integração na em- 
presa, 

Considera-se, poiz, necessário e 
oportuno lançar uma acção que de 
algum modo possa contribuir para 
tal objectivo, pondo-se à disposição 
das PME's, praticamente sem so- 
brecarga para os respectivos orça- 
mentos, durante o período do es- 
tágio, técnicos com importante 
formação de base, capazes de de- 
sempenharem missões de apoio à 
gestão daquelas empresas, 

O que, por outro lado, permitirá 
também aprofundar o conhecimen- 
to das situações de desemprego e 
subemprego dos quadros técnicos 
e das carências das PME's neste 
domínio. 

Deste modo, o Instituto de Apoio 
às Pequenas e Médias Empresas 
Industriais (TAPMET), como orga- 
nismo especializado no apoio e 
assistência às PME's industriais, e 
os servicos competentes da Secre- 
taria de Estado da População e 
Emprego (SEPE), especializados 
em matéria de política de emprego 
e formação profissional, desenvol. 
verão as acções necessárias por 
forma a dar a melhor execução a 
um plano experimental visando a 
formação e integração empresarial 
de quadros, que se designará por 
FIEQ. 

2 — Às empresas candidatas a 
este apoio deverão ser selecciona- 
das tendo em atenção as seguintes 
prioridades e de acordo com a re- 
partição regional adiante referida: 


1) PME's objecto de apoio es- 
trutural por parte do IAPMEI; 
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2) PME's beneficiárias de finan- 
ciamentos da SEPE ou de avales 
do TAPMEI com planos de recupe- 
ração; 


3) Outras PME's credenciadas 
com necessidade de admissão de 
quadros. 

9 — €) Serão apenas considera- 
das solicitações apresentadas por 
PME's, conforme indicado no n.º 2, 
relativamente aos titulares de li- 
cenciaturas ou bacharelatos minis- 
trados pelas diversas escolas ofi- 
ciais do País, com interesse nos 
domínios tecnológico, de gestão e 
organização e economia, 

b) Poderá também ser incluída 
no âmbito deste plano experimen- 
tal a realização de estágios para 
possuidores de cursos especializa- 
dos, nomeadamente de design, em- 
bora neste caso se deva adoptar 
um esquema de funcionamento 
adequado à sua especialidade, 
eventualmente com a colaboração 
das respectivas escolas no âmbito 
de acordos a celebrar. 

c) A selecção e proposta dos 
candidatos é da responsabilidade 
das PME's, não podendo ser admi- 
tidos, por empresa, ao abrigo deste 
esquema de apoio, mais do que 
dois estagiários de entre as áreas 
referidas na alinea a). 

4 -— a) Este plano experimental 
compreenderá setenta e dois está- 
Eios igualmente distribuídos pelas 
Regiões Norte (Porto), Centro 
(Coimbra) e Sul (Lisboa). 

b) Para os estágios previstos no 
n.º 3, é fixado o número máximo 
de dezoito, no âmbito deste plano 
experimental; 

c) A selecção dos estagiários pe- 
las PME's deverá circunsecrever-se 
a recêm-licenciados em situação de 
desemprego pós-escolar comprova- 
da e ter em conta, entre outros, os 
seguintes critérios: 


1) Maior ambiguidade de forma- 
tura; 

2) Maior nota de curso; 

3d) Maiores encargos familiares. 


5 — À coordenação do FIEQ ca- 
berã a uma comissão coordenado- 
ra, funcionando junto do IAPMEI, 
a qual será por dois re- 
presentantes, a temmo completo, 
dos serviços da SEPE, e outros 
dois do IAPMEI, cabendo a presi- 
dência a um dos representantes do 
Instituto. 

6 -— O estágio decorrerá na em- 
presa, terá a duração de seis me- 
ses e incluirá, além do trabalho a 
realizar na mesma, dois cursos de 
formação com a duração máxima 
de três semanas, sendo o primeiro 
subordinado ao tema genérico 
«Empresa: técnicas gerais de ges- 
tão». 


T— a) O acompanhamento dos 
estagiários será feito, na medida 
do possível, durante a sua perma- 
nência nas empresas, pelas equi- 
pas de assistência do IAPMEI em 
colaboração com os quadros e di- 
rigentes daquelas. 


b) Os estagiários deverão apre- 
sentar um trabalho à sua escolha, 
sob orientação da comissão coor- 
denadora, relativo a um problema 
epecífico da empresa, 


8 — a) Cada PME beneficiária 
do FIEQ deverá celebrar com os 
estagiários um contrato individual 
de trabalho, obrigando-se à sua 
inscrição na Previdência Social e 
ao seu seguro por acidentes de 
trabalho, encargos estes que serão 
de sua conta, 

bj A cada PME será atribuido, 
durante o período de estágio, e por 
estagiário, um subsídio fixo de 
10 000500 mensais, não podendo, 
em caso algum, a remuneração do 
estagiário ser inferior a este mon- 
tante, 

Os encargos respectivos serão 
suportados pelo IAPMEI (40 %) e 
pelo Fundo de Desenvolvimento da 
Mão-de-Obra (FDMO) (60 56), fi- 
cando a cargo do Instituto o apoio 
administrativo mecessário. 

9 — O IAPMEI assegurará o pa- 
gamento da totalidade dos encar- 
gos para com terceiros, após o 
visto da comissão coordenadora, 
transferindo o FDMO para o Ins- 
tituto a parte da verba orçamen- 
tada que lhe corresponda suportar. 
No final da execução do plano ex- 
perimental será efectuado encon- 
tro de contas mediante apresenta- 
cão pelo IAPMEI da respectiva 
documentação. 

10 — Todos os acordos e contra- 
tos com terceiros serão celebrados 
por intermédio do IAPMEI. 


11 — Serã elaborado um relató- 
rio final pela comissão coordena- 
dora, do qual deverá constar uma 
análise geral da situação do de- 
semprego e subemprego dos qua- 
dros técnicos, a avaliação dos re- 
sultados deste plano, um mapa 
com as despesas efectuadas e pro- 
postas de acção subsequentes de 
âmbito mais geral, enquadrando 
iniciativas desta natureza ma polf- 
tica de manutenção e criação de 
postos de trabalho e formação pro- 
fissional. 

12 — Os membros da comissão 
coordenadora, a nomear no prazo 
de quinze dias pelo TAPMET e pe- 
los serviços competentes da SEPE, 
deverão comecar por elaborar uma 
programação detalhada da forma 
de execução e de acompanhamento 
deste plano. 


(Continua na pdg. 201) 
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Prática de Krigagem 
FERNANDO HUMBERTO MUGE (*) 
RESUMO 


Introduz-se a noção de estimador linear optimizante — estimador de Krigagem. Constroi-se o sistema 
de Krigagem para o caso estacionário, Dão-se exemplos de variáveis regionalizadas em que não é válida 
a hipótese de estacionaridade e definem-se para estes casos os estimadores de Krigagem Universal. 

Dão-se exemplos de Krigagem de variáveis intercorrelacionadas (Cokrigagem) e constroi-se o res- 
pectivo sistema de Cokrigagem. 

Indica-se a estrutura búsica de um programa de Krigagem. 


ABSTRACT 


It is introduced the notion of best linear unbiased estimator — Kriging estimator. It is presented the 
Kriging system of equations for the stationary case. 
They are given examples of regionalized variables for the non stationary case and are defined, for 
these cases, the best linear unbiased estimador — Universal Kriging estimator. 
- They are also given examples in which are kriged intercorrelationated variables (cokriging) and is 
presented the respective cokriging system. 
Ft is shoum the basic structure of a Kriging program. 
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1 — POSIÇÃO DO PROBLEMA 


O problema da Krigagem consiste em encontrar o melhor estimador linear possível de uma caracte- 
rística desconhecida a partir da informação disponível. 
A informação disponivel consiste geralmente numa informação experimental, por exemplo nos teores 


(*) Assistente do Instituto Superior Técnico; Investigador do I.N.L C no CV. R.M.U. T.L. 
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de troços de sondagem implantados no interior e no exterior do painel a estimar e numa informação estru- 
tural, por exemplo o variograma da Função Aleatória representativa da regionalização dos teores estudados. 

Sendo Z (x) a Função Aleatória estudada, pontual ou considerada como tal à escala da estimação 
projectada, estacionária de 2.º ordem, de esperança: 


E[Z(x)]) = m (constante) 


Vamos supor conhecido o seu momento de 2.º ordem: 


E ([Z(x)-m] [Z (x+h)-m]J=C(h) — covariância centrada 
ou E([Z(x+h)-Z(x)]2) =2. y(h) — variograma 
— A informação experimental é constituída por um conjunto discreto de N teores Z, (com i=1,N). 


E 


1 
Estes teores Z, são, ou teores pontuais Z(x,) ou teores médios Z, = | Z(x) dx definidos em suportes 


Cid Yi 


de dimensão v, (volume, superfície ou comprimento). Introduzindo a hipótese de estacionaridade todos estes 
suportes têm por esperança m: E|Z,]=m Y.. 
— O propósito da Krigagem é o de obter o melhor estimador linear do valor médio de Z(x) num domi- 


nto V (bloco, painel, totalidade do depósito ou simplesmente num ponto de espaço): Z, — £ (x) dx 


l 
V v 
“ 
qe “xr 
O estimador de Krigagem: Z, tomará então a forma: 4, = > RA 
l=1 


em que os X, são os ponderadores a afectar à informação experimental. 


2 — EQUAÇÕES DE KRIGAGEM 


O estimador de Krigagem deverá ter as seguintes propriedades: 
sã deve ser um estimador não enviezado isto é: E [Z a pe = 0==> > ki = 9 
| E 
| 


— A variância de estimação (variância de Krigagem) deverá ser mínima, o que conduz ao seguinte 
sistema linear de equações (sistema de Krigagem): 


Preta, myj-g= Tt V) AR DA UM 
J 


2) 


É 


o 


A variância de Krigagem escreve-se; 


ci =E[Z,—-Z;])'=C(,W+E— pa dk E (nm, V) 
] 


Podem-se estabelecer as equações de Krigagêm partindo do variograma y(h) em vez da covariância 
C(h). Obtém-se então o seguinte sistema de equações: 


Dn, n)FE= (vv) co Wi=a WB 
j 


>, 4 — 1 


nm, WD+E—7 (VV) (1) 


1] 


q 
| 
-[ 
| 
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Em notação matricial o sistema de Krigagem escreve-se: 


(6) x [3) = [4] —> 


7 (iva) eso 7 (Vi, Yi) Pd 7 (Vas VN) 1 
a] : 7 (Vi, Vi) Dela Cras 7 (Yi) reina: y (mi, Vw) 1 
ga, VD essas y Ca, cce y(Vy Vx) 1 


Note-se que a matriz de Krigagem é simétrica: 


A variância de Krigagem o; não é senão o produto da matriz transposta [A]* pela matriz 2.º membro 
[M] menos o termo Y (V, V). 


OBSERVAÇÕES AO SISTEMA DE KRIGAGEM 

1-— O sistema de Krigagem é um sistema de (N+1) equações lineares a (N+1) incógnitas que são 
os N ponderadores à das N informações experimentais Z e o parâmetro de Lagrange uy. Demonstra-se 
que o sistema de Krigagem é sempre regular, isto é, fornece sempre uma solução única. 

2— A Krigagem, estimador linear e não enviezado, é também um interpolador exacto isto é, se o 
suporte V a estimar se confundir com o suporte Ye de uma qualquer das informações experimentais dispo- 
níveis, o sistema de Krigagem fornece: 


— Um estimador Z idêntico ao teor conhecido Z, desse suporte v =V. 
— Uma variância de Krigagem nula: 9/= E( [Z, Ee Es b="0 


Isto significa que a «superfície Krigada» passa pelos pontos experimentais, 
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Esta propriedade não é verificada por todos os processos de estimação, em particular pelos processos 
de ajustamento polinomial utilizando o método dos mínimos quadrados. 


3— À escrita dos sistemas e variâncias de Krigagem com as notações simbólicas y ou C é perfeita- 
mente geral quaisquer que sejam os suportes Vo V, da informação e o suporte V da grandeza a estimar. 

Assim, certas informações v, podem estar inclusas em V: v.cV. As informações podem ainda sobre 
por-se parcialmente v. V; É o Pode-se ainda considerar como (N+1)* informação a média aritmética das 
N informações já consideradas. 

O tipo de estrutura subjacente caracterizada pelo variograma Yy(h), embricado ou não, isótropo ou 
anisótropo, também não interfere na forma genérica do sistema de Krigagem. 


4 —O sistema de Krigagem e a variância de Krigagem dão, além disso, conta dos seguintes factos 
essenciais que condicionam toda a estimação: 


— Geometria do domínio a estimar: termo Y (V, V) de o2. 
— Distância entre o estimado e o estimante: termos y(v,, V) do 2.º membro [MJ]. 


— Geometria interna da informação: termos y(v,, v,) da matriz do 1.º membro [G]. 
— Estrutura subjacente caracterizada pela função estrutural semi-variograma y(h). 


Note-se ainda um facto extremamente importante, a variância de Krigagem: 

e, = [lt. [IM] — CV, V) não depende senão da estrutura y(h) e das geometrias e não dos valores 
efectivos da informação. 

Pode-se portanto, conhecendo y(h), prever a qualidade da Krigagem deste ou daquele painel V com 
esta ou aquela configuração da informação. 

Este facto permite optimizar o reconhecimento das jazidas (V.p. secção Geostatística como instru- 
mento de Flaneamento Mineiro). 


o— A matriz do 1.º membro [G] não depende senão das geometrias relativas internas (vo, Vi) das 
informações disponíveis e não depende da geometria do painel V a estimar. Isto deve-se ao facto de o sis- 
tema de Krigagem ser um sistema linear de equações o que conduz a que a duas configurações de informa- 
cão idênticas corresponda a mesma matriz [G]. 

As duas matrizes-coluna solução [%] e [)'] obtém-se então pelo produto da mesma matriz inversa 


[G] dá pelas 2 matrizes coluna do 2.º membro: 
= —1 
[A] = [6] . IM] e [NX] = [6] . [M'] 


2.1. — Composição de Krigagens 


Na prática o estimador global de uma zona ou mesmo de todo o jazigo é obtido por composição dos 
estimadores de unidades locais (painéis ou blocos). 

Como vimos atrás (observação 5 ao sistema de Krigagem) o suporte V a estimar só intervém no 
2.º membro do sistema. Assim a solução do sistema de Krigagem: [A] = [G] a [M] é linear em V., 

Seja à, (x) a solução da Krigagem do teor pontual Z(x) implantado em x, a partir de uma dada confi- 
guração de informação (C). A. partir da mesma configuração (C), procura-se agora Krigar o teor médio: 


V — E v f (x) d x 


definido no suporte V. A solução À (V) desta Krigagem será então: 
de (V) == === A (x) d x » VIGO 
VW 


Há portanto combinação linear das Krigagens pontuais para constituir a Krigagem global de V. 
Donde se conclui que é indiferente Krigar directamente a média global para todo o jazigo ou proceder 
por composição de Krigagens locais. 
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Esta propriedade é extremamente útil pois permite calcular imediatamente, por simples ponderação, 
o melhor estimador do teor global de um jazigo a partir dos teores Krigados para cada painel em que ele 
foi dividido. 

A estimação painel a painel é deste modo perfeitamente equivalente à estimação global. 


Observações: 


— A variância de Krigagem: o? | = E | [Zip ms Zy) | não se obtém por simples combinação linear das 
variâncias de Krigagem pontuais. 

— O teorema de composição de Krigagens supõe que todos os estimadores são construídos a partir da 
mesma configuração de dados. 


3 — KRIGAGEM EM PRESENÇA DE DERIVA — KRIGAGEM UNIVERSAL 
3.1. — Introdução 


Considerâmos até agora variáveis regionalizadas Z(x) tais que E[Z(x)]=m= constante. 

Existem no entanto muitos fenómenos para os quais a hipótese de estacionaridade não é válida. 

Assim, por exemplo, a plataforma continental afunda-se progressivamente e em média um ponto mais 
afastado da costa apresentará uma profundidade maior. Esta deriva m(x) é além disso, o factor preponde- 
rante de variabilidade entre as V.K. [1]. 

Em certas zonas de um jazigo pode-se detectar a existência de uma deriva (empobrecimento progres- 
sivo nas zonas limítrofes do jazigo, por exemplo). E portanto necessário de tomar em linha de conta a 
deriva quando se pretende efectuar a estimação por Krigagem dessas zonas do jazigo. 


j “a 


sondagens disponiveis 


Assim, por exemplo, consideremos uma lentícula mineralizada cujos bordos terminam em bisel, Pode 
acontecer que se tenha que estimar a possança média de um painel de borda a partir de sondagens situadas 
fora do painel. 

Uma Krigagem estacionária como qualquer outro processo de estimação estacionário sobreavaliará 
inevitavelmente a possança média real do painel. É portanto necessário ter presente a deriva das possanças 
na direcção q de empobrecimento da lentícula. 

Em pgravimetria é usual distinguir duas componentes do campo gravítico: uma componente regional 
e uma componente local. 

Esta separação é utilizada nos problemas de interpolação e extrapolação da gravidade distinguindo-se 
uma componente determinista denominada deriva regional («regional trend») e uma componente não deter- 
minista que tome em linha de conta os desvios locais tomados em relação à deriva. 

Os métodos correntemente utilizados para resolver este problema são métodos de ajustamento polino- 
mial que consistem em minimizar a soma dos quadrados das distâncias entre os dados experimentais e a 
deriva (mínimos quadrados), O modelo subjacente a tais métodos pressupõe que o fenómeno é a soma de 
uma componente determinística e um termo de erroc(x) com média nula, variância constante e sem auto- 
«correlação (ruído puro): 


k 
2(x) = ba a, € (x) + elx) 
1=1 


Este modelo que pode ser adequado quando as variáveis aleatórias são independentes, revela-se con- 
tudo inapto para representar fenómenos espaciais em que a autocorrelação desempenha um importante papel. 
O modelo adaptado na Krigagem Universal (K, U.) é o seguinte: o fenómeno em estudo é visto como a reali- 
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zação de uma função aleatória Z(x) que é a soma de dois termos: Z(x)=m(x)+Y(x) em que m(x) é à 
componente determinista chamada deriva que é igual ao valor expectável de Z(x) no ponto x: 


m(x) = E[Z(x)] 
Admite-se que a deriva admite uma representação sob a forma: 
m(x) = > a, f! (x) 
l 


Y(x) é a componente aleatória denominada resíduo. Este termo não deve ser confundido com o termo 
c(x) dos mínimos quadrados e isto porque Y(x) é uma função aleatória com estrutura própria. 


3.2. — Sistema de Krigagem Universal 


1 
Suponhamos que queremos estimar o teor médio: Z, = = | v Z(x) dx a partir da informação dispo- 
1 . 
nível constituída, por exemplo, pelos N teores médios Z, = — | Z(x) dx em presença de: 
v, v 
i “é | 
k 
— uma deriva m(x) = > a, f(x) de forma suposta conhecida 
I=o 


— um variograma y(h) 
* 
Tomemos um estimador linear da forma: Z,., = > À; Z, 
Í 


Os ponderadores são escolhidos de modo a que as condições de não enviezamento e de variância mínima: 
sejam respeitadas: 


1) E[4, —Zyy 


] = Ú não enviezamento 


O que conduz a: 


1 1 
aa. -=— | f =0 
[> » Í. (x) dx Tl, o 4x | D a, 


O que é uma identidade em a,. Donde; designando por: 


| l | 
A | T" (x) Cu 
v 


| 
| S ab=b v1=0% 
i k 


vem: 


(2) Variância de estimação minima 


A minimização de aa = E[2,-20,]? sujeita aos constrangimentos de não enviezamento obtém-se 
classicamente através do formalismo de Lagrange, e conduz ao sistema de Krigagem seguinte: 


as ge am 
à MTP DS ud=70W M=1N 
| 


e , 
à Md, =b, 4, =0,k 


k 
[e l — 
A=SNADED mdb, =7(V,V) 
Í |=0 


144 TÉCNICA 451/452 


O 


Observações: 


— O sistema de Krigagem Universal é um sistema de N+K+1 equações lineares a N+K-+-1 incógni- 
tas que são os N ponderadores de Krigagem à e os (K+1) parâmetros de Lagrange. 
Demonstra-se que este sistema é sempre regular, isto é que admite uma solução única, se e somente 


se as (K+1) funções f!(x) forem linearmente independentes sobre o conjunto das N informações, quer dizer 
se o sistema: 


a re E os 
3 e Ra ER o q=0, 4, 
I 
— O sistema de K. U. conduz ao sistema de Krigagem estacionária fazendo K = 0: m(x) = a0 = cons- 


tante. Com efeito f (x) = 1 e bo = 1, V, Dos (K+1) constrangimentos de não enviezamento: 


| 1 = 
>. A; b,, = by MW, só resta um: b> À, =, 
| | 


— O sistema de K. U. aparece como um sistema de Krigagem estacionária ao qual se juntaram K equa- 
ções suplementares correspondentes aos K termos a, f'(x), 10, da deriva m(x). 

A matriz 1.º membro [GU] do sistema universal obtém-se completando a matriz 1.º membro [G] do 
sistema estacionário com K linhas e K colunas: 


Ea — 
ss aaed 1 | K 
Viv, v,) AP E Y (vs Va) PRP E ee eia 1 De, edi CI UA CS DS 2 
[K] | - | 

1 1 ERRAR E: A 1 k 
by, Db, D,, 
| 7(vs E) ESC TR Qn ái y (o a? ad be, puro pen ua Ane b.. q add fia: Pol La by. 
[GU] = | a a a pes ad Desnsinorannçintea 0 
: “DR j PDS RD 0 

A sic pn pio 0 | o 
a A 0 Dao Rin cao 0 


O e T E ia ua 0 asd O asno 0 
vi “gm aa 


Note-se que a matriz [GU] permanece simétrica, 
Para as matrizes coluna desconhecidas [AU] e 2º membro [MU] vem também: 


A, 7 (vs V) t 

a 7 (v, V) N + 1 linhas 
RUJ=| [MU] =| 1 | 

x E k linhas 

A bh, | 
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O sistema de Krigagem Universal sendo: [GU]. [k U] = [MU], a variância de Krigagem escreve-se: 


sky = [AU]. MU) = (V,V) 


— Da Krigagem estacionária a esperança desconhecida à Krigagem em presença de deriva vão apare- 
cendo condições de não-enviezamento cada vez mais numerosas. A esperança desconhecida mas suposta 
constante m, (Krigagem estacionária) aparece a condição única: jT À, = 1. Quando a esperança é desco- 


| 
nhecida e não estacionária m(x), aparecem (K+1) condições sobre os ponderadores. 

Quanto mais numerosas são as condições, quer dizer quanto mais complexa é a estrutura da espe- 
rança desconhecida, mais difícil será a estimação, quer dizer mais forte será a variância de Krigagem. O que 
se ganha protegendo os riscos de enviezamento devidos a uma eventual deriva perde-se em variância de 
Krigagem., 

4. — COKRIGAGEM 
4.4. — Introdução 

O objectivo da cokrigagem é o de melhorar estimações locais tomando em consideração a informação 
adicional trazida por uma variável diferente da que se pretende estimar. 

Acontece por vezes que se estima uma grandeza, o teor em Pb dum painel, por exemplo, não só a 
partir da sua regionalização em Pb, insuficientemente reconhecida, mas também à custa doutra regionali- 
zação, a dos teores em Zn por exemplo, 

Um outro exemplo é a estimação do teor em U O, à custa da informação existente sobre teores e 
da radiometria dos furos de sondagem o que constitui outra variável regionalizada cuja estrutura apresenta 


forte correlação espacial com a dos teores em U O,. 
É evidente que a utilização da cokrigagem só é útil quando: 


— a informação disponível na variável em estudo é insuficiente para obter boas estimativas; 


— a correlação espacial entre as variáveis é suficientemente forte para que a variável adicional traga 
algum ganho de informação que sirva para a estimação da primeira. 


Sejam então K funções aleatórias intercorrelacionadas;: 
Z (x), EX), ..., Zp(x) 
sendo os seus momentos dados por: 
momento de 1.º ordem: E (Z,(x)) = mM, 
momento de 2." ordem: E(Z,, (x+h) Z,(x)) - mom = Cy (h) 


OU: 
E([Z, (x+h)—-Z, (x)] [2 (x+h) —, (x)1) =2 Vide (h) 


Ya (h) designa-se por variograma cruzado das variáveis k com k. 


4.2. — Sistema de Cokrigagem 
Consideremos então a estimação do teor pontual Z, (x, ) a partir da informação constituida por: 


N, teores pontuais Z (x,), ik=1, N, A, 


* 
Seja Z, (X.,) O estimador de cokrigagem: 


R 
ol) = Da 20) + DD Ma Za) oe DM Ze Gi) 
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A condição de não enviezamento conduz a: 


B(ô, Cu) Ás Cu))=0 => Di=1e SD m= 
ik 


1k 


n 


A minimização da variância: 


or, E E( [2x (Ho) “A, (X%,) ] ) 


V k&k, 


com os K constrangimentos precedentes conduz ao sistema de D> N | +K ) equações lineares designado 


por sistema de cokrigagem: 


[ES 
A + dk Tae (io Ho) = Tok! (X Xv) Wk ex gk 
k=1 dk=1 
Ny, 
2 Mix, =1 
a 
No 
DM =0 4 k5k, 
jk=1 


A variância de krigagem é dada por: 


To — py > e Tete Co Ha) Tp — Trk Co %) 


k dik 


Supondo duas variáveis 1 e j vem em notação matricial: 
[CK] X [A] = [M] —> [DJ =[CKI] 


A matriz do 1.º membro ou matriz de cokrigagem escreve-se: 


Ti Cp E) ——— (sp exi) Rc aj! 
| 
Ti Ce K)———— YOtsçis Xwça) FT Otsgio Yy) 
Fã Yyp E) ———>+, (9 4p Xw;1) Y; (Yi Y1;) 
[CK] = 

Tiso Rj) ——— (Yes Ka) Os; Yyj) 
RR 0— 
O Li) 1 
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"x [MJ] 


Yi Oto Yxj1) 


a (Yyp YNj) ) 


a 3 (Yxjho Yx;1) 
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e as matrizes coluna solução [)] e matriz 2.º membro [M] são: 


di Ti (X,, V) 

Axl +, (Xi V) 
e ds; M = LARA V) 

hj) Yi tYsy v) 

Hy 1 

H O 


A variância de krigagem vem na mesma notação: 
cer = —Yy (VV) + Dt x (M] 


fórmula esta idêntica à de krigagem estacionária. 


5. — ESTRUTURA BASICA DE UM PROGRAMA DE KRIGAGEM 
Os dados de partida de um programa de krigagem são os seguintes: 


— um ficheiro de amostras com os respectivos teores e coordenadas. 
— um ficheiro de blocos a serem estimados (características geométricas e implantação espacial). 


— um modelo para a função variograma, 
Para cada bloco a krigar será repetido o processo seguinte (ver fig. 1): 


— Considera-se um bloco. 


— O ficheiro das amostras é percorrido, procurando-se as amostras que caem na zona de influência 
(auréola) do bloco. 


— Uma vez encontradas, são calculadas as covariâncias destas amostras, o que define a matriz 
1.º membro do sistema de krigagem [G). 


— São calculadas as covariâncias destas amostras com o bloco a krigar o que define a matriz coluna 
2.º membro do sistema de krigagem M (ver fig. 1). 


— O sistema de krigagem: [G].[)] = [M] é resolvido para a determinação da matriz-coluna dos 
ponderadores ou coeficientes de Krigagem [)J]. 


— O estimador de krigagem do teor do bloco é obtido multiplicando os ponderadores pelos teores das 
amostras retidas. 


— À precisão da estimação é calculada pela fórmula (1). 


5.1. — Observações 


5.1.1. — Pesquisa das Amostras 


Na maior parte dos programas que utilizam os métodos clássicos (método dos polígonos, das distân- 
cias ponderadas, etc...) a pesquisa das amostras que intervém na estimação dum bloco é realizada a duas 
dimensões por polígonos ou por cones no espaço a 3 dimensões. Isto é feito para assegurar uma represen- 
tação uniforme para todas as direcções do espaço de modo a evitar o efeito de «écran» que se verifica 
quando ocorrem concentrações de amostras preferencialmente em certas direcções. 


Este cuidado é desnecessário em krigagem e isto porque a introdução das covariâncias entre as 
amostras provoca uma distribuição automática da influência pelos constituintes do grupo evitando a sua 
sobrerepresentação no sistema de krigagem. 
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PROGRAMA DE KRIGAGEM 
(DIAGRAMA DE BLOCOS) 


INPUT: | 
EGOR DENADAS DA INFOR 
MAÇÃO 


"SR AMETROS DO VARIO- 
GRAMA 


DIMENSOES E ORIENTA- 
RO DOBLOCO A KRIGAR 


CALCULO DO VALOR 
MEDIO DO VARIOGRA- 
MA NO BLOCO (VARI- 
ÂNCIA DUM PONTO 

NO BLOCO) 


PESQUISA DO BLOCO N a FiM DO 
A KRIGAR | cido 


PROGRAMA 


PESQUISA DA INFOR- 
MAÇÃO EXISTENTE 
NA AURÉOLA DO 

BLOCO 


CALCULO DAS COVARI-| 
ÂNCIAS ENTRE AS 

INFORMAÇÕES (EX: 
SONDAGEM A SONDAGEM 


CALCULO DAS COVARI- 

| ANCIAS ENTRE AS IN- 
FORMAÇÕES E O BLO- 
CO A KRIGAR 


RESOLUÇÃO DO SISTE- 
MA DE MAT HERON 


DEFINIÇÃO DO ESTIMA- 
DOR DO TEOR DO BLOCO 


Fig. 1 


5.1.2. —- Cálculo das Covariâncias 


A covariância do bloco Z com a amostra X.: q, « é por definição: 
E DR 


=C+C- V | | | Y (v (x—x.)2 + (=) + (2—2,)? | dx. dy. dz 


em que V é o volume do bloco Z ex, y, e z são as coordenadas da amostra X.. 


O cálculo desta covariância é feito pelo processo numérico ilustrado na fig. 2. 
Os valores médios y(v,, v,) são calculados numericamente, discretizando cada um dos suportes v,, V,, 


o que conduz a assimilar cada suporte v, à reunião de N, pontos x, k = 1, N, dispostos segundo malha 
regular em V. 


Fig. 2 
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Uma nova parametrização de jazidas minerais 


RESUMO 

Apontam-se algumas limitações à prática habi- 
tual do critério de rejeição baseado sobre pequenos 
blocos e respectivo teor tal-qual. 


É-se assim condusido à consideração do bloco de 
beneficiamento e à construção da característica opti- 
mal, a qual permite orientar a escolha do projecto 
optimal e, traduz, portanto a parametrização mais 
significativa. 


Mostra-se seguidamente que o controlo do des- 
monte é realizado após parques de homogenização e 
desenvolve-se a teoria básica que permite a avaliação 
das políticas de constituição e retoma de tais parques. 


1. INTRODUÇÃO 


1.1. Só em fase preliminar do planeamento da 
exploração de uma jazida se justifica reportarmo-nos 
à distribuição de teores exibida pelos troços dos tes- 
temunhos das sondagens S, visto que o respectivo 
histograma tem maior variância de dispersão do que 
o dos blocos tecnológicos V, sobre os quais se irá 
exercer a decisão de tratar ou não tratar ou ainda 
de como tratar o minério, 

A conhecida relação de aditividade de Krige 
mostra como se relacionam as variâncias de disper- 
São em causa (1). 

ke) z | 
von “anT ay 

Deve também notar-se que se sabe calcular “A 
uma vez que, pelo menos, exista variograma e V 
seja uma partição de D ou ainda V <£D: 


CT PE e Does. Mega 
vn = UD,D) — NV, V) 


Por outro lado, se as sondagens forem suficien- 
temente abundantes, os primeiros momentos espa- 
ciais das duas distribuições serão aproximadamente 
iguais: 


JOSÉ QUINTINO ROGADO 
Universidade Técnica de Lisboa 
Instituto Superior Técnico 


ABSTRACT 

Some limitations are pointed out concerning the 
usual cut off criterion based on small blocks and 
respective crude ore grades, 

By this via the concepts of beneficiation block 
and optimal characteristic are built up. This charac- 
terístic is useful for framing the design of the opH- 
mal mining project. Thus, it às the convenient 1007 
for the parametric description of an ore body. 

It is also shown that the exploitation control is 
usually carried out on the output of a stockyard. 
The theoretical basis which governs the adequate 
stacking and reclaiming stockyard policies às deve- 
loped. 


Z, = Z, 
Ou, consequentemente, admitida a ergodicidade 
do primeiro momento: 


E(Z,] = BZ, 


Todavia, a estabilidade da lei de distribuição não 
se pode garantir (2), excepto no caso gaussiano que 
é anormal entre as variáveis regionalizadas, como se 
sabe. 

Quer dizer que o histograma dos teores verda- 
deiro dos blocos V é, na maioria dos casos, inaces- 
sível, 

O que se conhece é a distribuição dos seus esti- 
madores, que, no que se segue, consideraremos esti- 

* 
madores de krigagem Z.. 


Assim, por definição de estimador não enviezado, 
como é o de kKrigagem, será: 


(1) Aliás, a esta mesma conclusão se pode chegar sem recorrer à citada relação de Krige [5]. 


(2) É costume admitir-se a estabilidade da lei ge: 
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4 E] | 
Sr) Dy < 02 
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Ea 
Contudo, este estimador sub-estima a variância 


de dispersão dos teores reais, como mostra a relação 
de atenuação na sua expressão prática: 


se sa SE 
DE (Di) + 
D,=E[(Z2,—m)'] 


(D,)'=E [( Zi —m) | 


com. 


Ê 


N 
q fo sã] E 
ey = El (Z, =) | a > “ky 

i 

Portanto, o teor de corte baseado nos estimado- 
res de krigagem dos teores dos blocos V pode, já por 
esta razão, não ter muito significado, 

Porém, durante a exploração, admitindo mesmo 
que o bloco de decisão continua a ser o que foi como 
tal adoptado no planeamento, surge informação adi- 
cional que resulta da amostragem de controle que 
precede o arranque nas frentes. 

Daí que aplicando mais uma vez a relação de 
aditividade, será, já que nova krigagem é possível: 


Dy=(D,) + oxw 


e porque é certamente 


ter-se-á: 


(Diet (DS 


Isto é, o teor de corte inicialmente escolhido tem 
pouco interesse ao nível da exploração, porque ha- 
verá que corrigi-lo à medida que a informação nova 
surja, o que não tem significado, já que o critério 
de rejeição é um conceito global e não local, 

Para sair deste impasse, pode recorrer-se à kri- 
gagem disjuntiva com a qual, em situações favoráveis, 
é possível prever o histograma futuro dos teores 
dos pequenos blocos, condicionado pela informação de 
que se dispõe ao planear a exploração, 

Não abordaremos aqui este importante progresso 
da Geoestatística mineira. 

1.2. Por outro lado, o teor de corte baseado no 
teor tal-qual é, em geral, uma aproximação grosseira 
da análise económica que assume, posto que a mar- 
gem monetária por unidade de peso de minério tal- 
-qual calculada à saida do concentrador é: 


by=Rypey— (m+Dy 


= nypZy — (m +, 


—— —=—— ———————— — 


Daí, que o teor de corte económico (baseado no 
critério marginalista mineiro) se não possa estabe- 
lecer com o teor tal-qual, embora razões tecnológi- 
cas, muitas vezes, exijam um primeiro nível de rejei- 
ção baseado naquele teor. Com efeito, a maioria dos 
processos de tratamento não pode garantir teor mé- 
dio de concentrado acima de certo limite se o con- 
centrador não for alimentado com teor médio tal- 
-qual superior a determinado limiar. 

Surge, todavia, desta análise, uma situação nova 
que importa explicitar. 

Efectivamente, o rendimento em metal mw O 
rendimento em peso R e o teor do concentrado c, 
determinados por ensaios padrão bem definidos (1), 
quando reportados a resultados industriais, são variá- 
veis regionalizadas mais regularizadas do que o teor 
tal-qual. 

Quer dizer, os variogramas respectivos têm am- 
plitudes superiores âquela que exibe o de Z.. 

E por essa razão, induzem volumes de influência 
maiores. É, portanto, mais fácil encontrar similitudes 
utilizando parâmetros mineralúrgicos do que quando 
se lida com teores tal-qual para agrupar blocos for- 
mando «clusters» |]JV coerentes que constituirão os 
blocos de beneficiamento, 


Contudo, os processos industriais de que os refe- 
ridos ensaios padrão pretendem ser as réplicas não 
são, como se sabe, operações totalmente determinís- 
ticas. Contêm, por isso, componentes aleatórias, cer- 
tamente correlacionadas, em certa medida, com Z. 

Quer isto dizer que os variogramas das variáveis 
reais m, R, e c, devem exibir significativo efeito 
de pepita. 

Portanto, se se pretender estimar mn, R ec, 
tomando por suporte o volume V, que é normalmente 
pequeno, pois se trata do bloco de escolha de des- 
monte, é de esperar que as variâncias de krigagem 
aumentem sensivelmente em relação às dos teores 
tal-qual e, além disso, que o chamado efeito de 
«ecran» se desvaneca acentuadamente, tornando os 
cálculos mais onerosos. 

Deve, por isso, fazer-se a estimação baseada em 
blocos de volume [JV nitidamente superior aos usados 
quando se estima os teores tal-qual. 

Todavia, a avaliação da massa mineral e do res- 
pectivo teor médio tal-qual, é óbvio que deve fazer-se 
à custa dos dados directamente obtidos na prospec- 
cão, pesquisa ou reconhecimento. 

É de notar, porém, que muitas vezes não se 
ganha grande coisa em krigar, sendo preferível efec- 
tuar os cálculos baseando-nos nos volumes de influên- 
cia: a variância da extensão à jazida do teor médio 
das sondagens conduzirá a erro razoável para o teor 
médio e massa mineral assim calculados. 

1.3. Outro aspecto que é relevante, embora pou- 
co referido, é que o histograma da jazida ou de parte 
importante da mesma não é estável às operações de 
desmonte, se se não tomarem cuidados apropriados 


(1) Para cada tipo de minério que se haja caracterizado na jazida. 
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Daí, a vantagem, ou mesmo a necessidade de 
organizar a exploração tentando estabilizar o pri- 
meiro momento do teor do concentrado e, se possível, 
também o do rendimento em peso ou em metal res- 
pectivos, o que se consegue se se puder explorar ré- 
plicas da jazida ou de parte importante desta refe- 
ridas a certas bases de tempo [2]. 


2. CARACTERÍSTICA OPTIMAL 


2.1. Entende-se por parametrização o estabele- 
cimento de bijecções entre a tonelagem tal-qual ou a 
massa mineral ou ainda o teor médio e determinado 
parâmetro de contrôle que pode ser o teor de corte 
do minério tal-qual, o teor de corte do concentrado, 
o valor de corte do rendimento em metal, etc. 

Quando a parametrização se faz à custa das 
variáveis mineralúrgicas, põe-se o problema de esco- 
lher, para cada bloco de beneficiamento |JV, o mais 
adequado ponto de funcionamento de forma a maxi- 
mizar a margem monetária total não actualizada EB, 
minimizando a perda em metal, uma vez fixadas a 
reserva total a explorar Q e a taxa de desmonte r. 

Tem-se então: 


B=Up-(fC+M(Q,r) + L(Q,r)) 


Considere-se, para resolver este problema, o es- 
paço das configurações mineralúrgicas TX U. 

Trata-se, obviamente, do espaço dos estados do 
sistema jazida/beneficiamento, 

Cada afixo Pc TX U representa um ponto pos- 
sivel de funcionamento desse sistema, que resultará 
de explorar todos ou parte dos blocos de beneficia- 
mento |V, realizando, em cada um deles, determi- 
nada operação mineralúrgica (considera-se, por ra- 
zões metodológicas, como operações distintas, regu- 
lações diferentes do mesmo processo mineralúrgico). 

Notar, que este espaço não inclui componente 
temporal, embora se possa falar em trajectória do 
afixo P, significando com isso que se adopta certa 
relação de ordem nas configurações possíveis do sis- 
tema. 

Esta é, aliás, uma situação típica da programa- 
ção dinâmica, quando aplicada à estrutura de um 
agistema susceptível de combinatória entre estados. 

Assim, o estado inicial do sistema será definido 
pelo afixo P(T,U,): O estado inicial é, portanto, 
a imagem do minério in situ. 

Muito raramente, como á sabido, esta configu- 
ração inicial é aceite pela tecnologia de juzante, 
quanto mais não seja porque T é inferior ao mínimo 
por ela exigido. 

Há pois que considerar outro afixo P(T,U) para 
que aquele pre-requisito seja satisfeito. 

Por isso se põe a questão de saber como organi- 
zar a exploração do sistema jazida/beneficiamento, 
de forma a passar de P a P com perda mínima em 
metal U, uma vez que se não conhecem processos de 
exploração sem perdas. 


E óbvio, que a perda em metal será minimizada, 
se de P se passar a P seguindo a trajectória extre- 
mante (maximizante) da funcional linear 


=[lsilar 


Trata-se claramente de um problema do cálculo 
das variações e, como tal, pode ser resolvido no do- 
mínio do continuo, 

Usar tal método implica, porém, interpolações 
analíticas das aplicações R — c trabalhosas e, como 
sempre, duvidosas, sem vantagem adicional, por- 
quanto a classe dos processos mineralúrgicos prati- 
camente admissíveis forma conjunto discreto finito. 

Assim sendo, poder-se-ia, entre outros, usar O 
algoritmo de Bellman discreto, pois é válido o prin- 
cípio da optimalidade. 

Felizmente porque R(7) é um convexo, mostra- 
-se [3] que o conjunto de operadores 


dr 
“LARS 


construidos transferindo sucessivamente cada bloco 
de beneficiamento UV do ponto de funcionamento 
correspondente a concentrado de teor mais baixo 
(e, portanto, de rendimento em peso mais elevado) 
para o ponto de funcionamento possível relativo ao 
teor imediatamente a seguir, fornece critério muito 
mais simples e muito menos oneroso para maximi- 
zar 1. 

Com efeito, estabelecida uma relação de ordem 
total em E, induz-se uma ordenação em (T). 

Assim, fixado T, basta procurar a PD (E) CE 
ordenada, cujo supremo conduz a atingir-se T. 

Demonstra-se [3] que, deste modo, se alcança 'T 
com perda em metal U minima: 


* 
A R 7] 
* (Re) 


Variando T, obtêm-se a bijecção optimal de U 
em C, a qual é, aliás, o convexo superior das confi- 
gurações possíveis do sistema, projectado no espaço 
UXC. 

Aquela bijecção chama-se característica opti- 
mal [4], porque traduz a curva de Lasky generali- 
zada optimal da jazida, tomando em consideração o 
desenvolvimento tecnológico coevo e a informação 
disponível, 


nas [map PA 


2.2. Em relação a esta característica optimal 
da jazida, qualquer projecto de exploração fica de- 
finido pelos parâmetros 


M(G,r) +L(Q,r), que determinam a taxa de 
desmonte r em minério tal-qual a Q 
fixado, ou vice-versa; 


(1) Notar que, no espaço U XC £ representa uma recta de coeficiente angular f/p e de traço em U igual a (M+Q)/p. 
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p que reflecte o mercado; e 
f que introduz a localização da jazida atra- 


vês o frete. 


A margem monetária global não actualizada de 
qualquer projecto é expressa, no espaço U X €C, por 


B/p = U-Éf 
1 
EC =0t/p+(M(Qr) +L(Qr))/p 


onde a parcela subtrativa representa os custos, ex- 
pressos em unidades do valor pro-rata de metal con- 
tido no concentrado e B/p a margem global não 
actualizada expressa na mesma unidade (2). 


Em relação a determinado projecto, o teor PT 
optimal é o que corresponde a Sup (U- £) e ao ponto 
de funcionamento (0,6) optimal, 

Notar que pode dar-se o caso de ser, por toda a 
parte, U- £ negativo, 

Num caso destes, não querendo funcionar com 
margem negativa (isto é, desenvolver actividade sub- 
vencionada), haverá que mudar de escalão de produ- 
ção, aumentando-o para diminuir M(Q,r) +L(Q, Fr), 
isto é, a despesa total de lavra e beneficiamento, 
fixado Q. 

Em certos casos, por razões de integração na 
macro-economia, poderá não ser possível ou conve- 
niente aumentar o escalão de produção. Haverá en- 
tão que rejeitar blocos da jazida, mantendo r fixo. 

Mostra-se [3] que a política que, neste caso, 
minimiza as perdas em metal, consiste em abandonar 
prioritariamente os blocos cujos teores dos concen- 
trados são os menores. 

Ao seguir-se tal política, Q deixa de estar fi- 
xado (porque diminui) e, portanto, mantendo-se r, 
M(Q,r) + L(Q,r) diminui aproximadamente na pro- 
porção da variação de Q. 

Nas aplicações, deseja-se conhecer M + L com 
o maior rigor possível, o que só se obtém elaborando 
ante-projecto ou projecto concreto. 

Pode então utilizar-se qualquer critério aceitá- 
vel de optimização da exploração, com Q livre, em- 
bora limitado superiormente à reserva avaliada mas 
com r fixado ou vice-versa, «mutatis mutandi». 

Encontrar-se-á assim (U- £) coerente com o 
critério adoptado. Contudo, se se quiser garantir, 
dentro de tal critério, perda mínima em metal, os 
blocos eventualmente a rejeitar deverão ser, priori- 
tariamente, os de menor teor de concentrado no ponto 
de funcionamento correspondente à característica 
optimal teórica. 


3. CONTROLE DE DESMONTE 


Já se referiu que assegurar, para determinada 
base de tempo, que o minério que se extrai da jazida 
seja uma sua réplica, é condição necessária e sufi- 
ciente para manter estáveis os histogramas, 


E claro que esta prática implica mistura de des- 
monte, a qual, por si só, nada diz sobre a variabill- 
dade do minério tal-qual extraído, considerada que 
seja outra base de tempo, 

Trata-se de ponto que merece análise mais deta- 
lhada. 

A variabilidade, em relação a um valor médio de 
qualquer parâmetro de controlo, avalia-se, como se 
sabe, pela variância, a qual só tem significado se se 
fixar: 

— O suporte (ou incremento) 

— o volume do lote (base de tempo) 

— à localização, na jazida, dos múltiplos in- 
crementos que irão constituir o lote. 


A importância da localização resulta da acessi- 
bilidade na jazida não ser total e, por isso, os incre- 
mentos que constituem determinado lote a controlar 
não podem ser colhidos de forma uniformemente dis- 
tribuída em toda a massa daquela, mas apenas em 
determinados painéis sucessivamente seleccionados e 
preparados para desmonte em consequência de cons- 
trangimentos vários. 

Admita-se, pois, que se escolhem n painéis V, 


cujos teores médios krigados são m, e que deles 
se extraem volumes v, na proporção p,, cuja reunião 
constitui o lote, a dispersão do qual se deseja contro- 
lar. Sejam ainda os v, partições de V, ou v, <€V.. 

Se v, por hipótese, puder ocupar qualquer loca- 
lização no painel V,, os estimadores das respectivas 
dispersões de teor serão 


| Ds) | = ev + AVoV)—r(n— vi) 


Portanto, o estimador da dispersão do teor do 


(3) 
lote L, será 


[Déo | -D(m8| Dé v)]): | E (1 Sã n) 


Por outro lado, o estimador do teor médio de L 


será, obviamente, 
. e 
| m(1.)| = p: Mm; 


Suponha-se também que os incrementos são se- 
quencialmente depositados num parque, adequada- 
mente manuseado, isto é, actuando como um opera- 
dor linear continuo (filtro linear). 

Seja Z o parâmetro em contrôle (para concre- 
tizar, Z pode supor-se um teor). 

Ter-se-à um fluxo de entrada no parque que se 
caracteriza pela série cronológica Z (t) e um fluxo 
de saída, pela série Z, (t). 

As massas ou volumes destes dois fluxos são 
iguais, em valor médio, calculado sobre base tempo- 
ral adequada, embora possam apresentar taxas «ins- 
tantâneas» distintas, 


(1) Notar que, no espaço U XC, £ representa uma recta de coeficiente angular f/p e de traço em U igual a (M+Q)/p. 


(2) Ou em unidades monetárias usuais, 


(3) L é, habitualmente, o volume ou massa da pilha constituída no parque. 
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rosa E = 


Não há embalagem que 
dispense um fecho inteiramente 
seguro. Para esse fim, é que 
existe Tesapack. À fita Tesapack, 
auto-adesiva, 
resiste nas 
embalagens de 
cartão fortes e 
espessas e 

mantém a sua aderência 
e elasticidade sob todas 


Tesapack 
um sistema de segurança 
para todas as embalagens. 


as condições atmosféricas. 
Para todas as embalagens, 
quer sejam grandes ou 
pequenas, largas ou estreitas, 
baixas ou altas, O sistema 
Tesapack - fitas e máquinas - 
oferece sempre um máximo de 
segurança. 

Tesapack aumenta ainda a 
rentabilidade das linhas de 
embalagem, quer com 


À 


| q 
máquinas simples, quer com 
máquinas automáticas de 
grande cadência. Tesapack, 
uma gama completa de fitas 
auto-adesivas para embalagens, 
em todas as indústrias. 


Tesa, sempre a solução do seu caso 
PARA MAIS INFORMAÇÕES CONSULTE O DEPARTAMENTO TESA DA BEIERSDORF PORTUGUESA, LDA, 


Est. de Barcarena - Queluz de Baixo Tel. 956171 ou Delegação no Norte - R. Gonçalo Sampaio, 1684-178 - Porto Tel. 6891153/4/5 
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Por isso, e por ser mais significativo, é preferi- 
vel substituir a variável tempo t pela variável q, 
significando massa ou volume: 


Z (1) = Z, (9) 
Z(t)—>2Z (lg) 


Se o parque actua como um filtro, existe a se- 
guinte relação entre as duas séries cronológicas: 


So q 
Za) as Ls (q) Ey (q-7)d + ) Z!q)rug-)dr 
“0 qAg 
cade r, q-7) é a resposta do parque à impulsão 
de Dirac centrada em q a qual, escoado Ag, prati- 
camente se anula 
Isto é, 


ala) =Z (q) * Tr 
é a regularizada de Z, pela função de amostragem r. 


Seja R.(h) =r, * r, o covariograma transitivo 
associado a r (q), será [1]: 


140) = (ne, (1) ) | ida) ala efa jagãa 
Dq Da 


Dn 


K 


onde: y (h) e y.(h) são os variogramas de Z e Z, 


respectivamente, calculados no suporte Aq 


| r.(g) rtg + hjda, por 
a PAN: 


% 
Rb) =1,* 1,5 
definição 


K é o valor médio do variograma y., ponde- 
rado por r (q). r (q') no suporte Aq; trata-se 
de uma constante que se pode determinar 
pela condição fronteira y (0) =o, isto é, 
K=y( +to)* R( +o): E, pois, o efeito de 
pepita de y, no suporte Ag. 


4. LISTA DE SIMBOLOS NÃO DEFINIDOS NO TEXTO 


Por outro lado, a variância pontual do teor cor- 
respondente ao lote IL, será (1) 


Di(o/Lj= 7, (LL)=y(h)! 
= vdh)* Ra(h)' 


= + RO 


Se a resposta do parque, à impulsão de Dirac, 
for estável, será R(h) constante. 


Em qualquer caso, pode saber-se calcular ou 


(2) L 
medir R(h) ' 
Donde 
z pl 1, 
D.(0/L) = D (o/L) * R,(h) 


Todavia, o que de facto se controla, por amos- 
tragem, é a variância de uma toma de volume s. 
Será então: 


D;(ejL.) = D,(0/L) — D,toje) 


= pan) — (hn) 


o que se sabe calcular, visto poder-se conhecer o 
variograma pontual do lote L à saida do parque, 
feita a desconvolução respeitante ao suporte Ap. 


Vê-se assim, que a questão do controlo, se reduz 
à determinação do variograma da série cronológica 
de entrada Z (t). Daí o interesse na fase de projecto, 
em utilizar as técnicas da simulação condicional, 
condicionada pela acessibilidade. 


“vip — variância do teor tal-qual verdadeiro do bloco V na jazida D 
Cs/p — idem da sondagem S 

Usfv — idem da sondagem 8 no bloco V 

Y(D, D) — variograma médio da jazida D 


Y0V, V) —— > + do bloco V 

MV Pod MM P; — » » do bloco Vou v, respectivamente 

Es Z. — média espaciais de A. e A. respectivamente, ElZ 1, E[Z] -., eSperanças de Am AR etc., 
respectivamente 

Cy — variância média de krigagem do bloco V 

B. — dispersão (variância à priori) dos teores verdadeiros 


(1) Notar que, tratando-se de variância pontual, o efeito de pepita será nulo. 
(2) Recorrendo, por exemplo a um analógico, método adoptado no Laboratório de Preparação de Minérios do 1. 8. T. 
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(D$)*, [DÍ (v,)]* — dispersão dos estimadores de krigagem dos teores de V e v,, respectivamente 


m = E |Z,.] — esperança dos teores verdadeiros 

Cy — variância de krigagem do bloco V, 

(DA ) ++ — dispersão dos estimadores de «segunda» krigagem 

BR. — rendimento em peso do concentrado obtido do bloco V 
Ny — rendimento em metal do concentrado obtido do bloco V 
p — valor da unidade do metal 

Cy — teor do encentrado obtido do bloco V 

(n+ — custos de lavra e tratamento por unidade tal-qual do bloco V 
C — peso do concentrado global 

U — metal contido do concentrado global 

f — frete por unidade de peso do concentrado 

M(Q,r) + L(Q,r) despesa global de lavra e tratamento, fixadas Q e r. 

É à — teor do concentrado global 

R — rendimento em peso do concentrado genérico 

fe — teor do concentrado genérico 

mo — estimador do teor médio do bloco v, 


p — proporção do bloco v, na constituição do parque L, 
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Teoria e prática de amostragem de materiais a granel 


segundo o formalismo de P. Gy 


RESUMO 


Explicam-se os mecanismos geradores dos dife- 
rentes erros de amostragem, identificam-se as varid- 
veis de que dependem, estimam-se os seus dois 
primeiros momentos e determinam-se as «qualidades» 
que uma amostragem deve possuir e como se pode, 
na prática, assegurar-lhas. São analisados os pro- 
cessos fundamentais de amostragem — partilha e im- 
plantação de incrementos — através de dois modelos 
matemáticos: Modelo Probabilista e Modelo de Inte- 
gração, utilizando a Estatistica Clássica e Geoesta- 
tistica. 

Apresentam-se exemplos de aplicação a materiais 
a granel da indústria mineira, nomeadamente con- 
cepção e cálculo de instalações automáticas de amos- 
tragem. 

Tem, contudo, maior generalidade, sendo aplicável 
a materiais a granel de outras indústrias, nomeada- 
mente metaliirgica, química, farmacêutica, agricola 
e alimentar, do cimento, do vidro, da cerâmica e do 
papel, q resíduos urbanos e industriais e ainda a 
materiais contínuos, sólidos ou líquidos, como chapas, 
barras, lingotes e produtos metalúrgicos fundidos em 
deslocamento. 


INTRODUÇÃO 


A amostragem tem uma importância fundamen- 
tal quer técnica (controlo de qualidade) quer eco- 
nómica (2) e é incontestavelmente uma realidade 
complexa. 

Efectivamente, a amostragem, no sentido lato, 
divide-se num certo número de andares, a saber; 


* andares de amostragem propriamente dita; 

* andares de fragmentação; 

* andares auxiliares (homogeneização, secagem, 
manuseamento, etc.), 


terminando pela determinação dum parâmetro de 
qualidade (análise química ou granulométrica, teor 


JOÃO CARVALHO MARQUES (') 


ABSTRACT 


The paper explains all types of sampling errors 


- and estimates their first and second moments. It 


defines the «qualities» and the rules to respect im 
sampling procedure. 

The main sampling procedures are studied 
through two mathematical models: Probabilistic and 
Integration Models, developed by Classical Statistics 
and Geostatistics. 

Applications to bulk materials of mining industry 
namehy predesign of automatic sampling plants. 

The models con be applied to bulk materials of 
many other industries such as Metalurgical, Chemi- 
cal, Food, Agriculture, Concrete, Glass, Ceramics and 
Paper, to industrial and city wastes and even to 
continuos solid and liquid materials such as ingots, 
rods, bars, steel plates and fused metalurgical pro- 
ducts in contínuos mooving. 


de humidade, etc.) objectivo final de toda a opera- 
ção de amostragem. 

Do ponto de vista teórico, só os andares de amos- 
tragem propriamente dita provocam erros de amos- 
tragem. Na prática, todos os andares sem excepção, 
todos os manuseamentos a que se submete o lote e 
suas amostras sucessivas podem introduzir erros, pois 
hã sempre o risco de alteração do material a amostrar 
(quer por perda de fracções, quer por adição de 
materiais estranhos, quer por modificações de estru- 
tura, quer ainda por erros involuntários e até volun- 
tários). Como o nosso objectivo é procurar criar um 
modelo geral de amostragem, desprezaremos provi- 


(1) Investigador do C. V. BR. M. U. L. e Docente do I. 8. T. 
(2) A Produção Mundial anual de matérias-primas minerais sólidas representa cerca de 2% do Produto Mundial Bruto. 
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soriamente todos estes erros operatórios, embora 
nunca esqueçamos a sua existência. 

A amostragem, no sentido restrito, será anali- 
sada nos seus dois processos fundamentais — o pro- 
cesso de implantação de incrementos e o da parti- 
lha — aplicados respectivamente a lotes não mani- 
puláveis e manipuláveis. 

O processo de implantação de incrementos será 
ainda estudado através de dois modelos teóricos: 
«Modelo de Integração» e «Modelo Probabilista». 
Cada um destes modelos interpreta a realidade con- 
soante a óptica, contínua ou descontínua em que é 
encarado o lote, resultante da ambivalência da natu- 
reza da matéria ou, se quisermos, do maior ou menor 
recuo com que analisamos a realidade. 

Veremos pois que, em cada andar de amostra- 
gem — amostragem no sentido restrito — podem 
ocorrer diferentes tipos de erros que, no caso mais 
geral, constituem o erro total. 


De forma esquemática indicam-se os erros que . 


integram o modelo geral da amostragem e cuja sim- 
bologia é a seguinte: 


ET -- Erro total 

EO Erros operatórios 

EM, — Erros do modelo 

EM, — Erros de materialização 

EC -- Erro de «découpe» 

EP -—- Erro de toma 

EQ Erro de flutuação de qualidade (a(x)) 


ED Erro de flutuação do caudal (erro de 
ponderação ) 

EI, — Erro de integração do termo regionali- 
zado 


EI, — Erro de integração do termo local 

El, ——Erro de integração do termo periódico 
EF -— Erro fundamental 

ES -—-Erro de segregação e agrupamento 
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Os erros de ponderação e de integração dos ter- 
mos regionalizado e periódico são estudados através o 
modelo de integração. 

Os erros de materialização, fundamental e de 
segregação são estudados no modelo probabilista, se 
bem que o erro de integração do termo local também 
seja analisado à luz do modelo de integração. 


1. ANALISE GERAL DO PROBLEMA 


1.0 — O modelo geral aplica-se aos andares de amostra- 
gem propriamente dita, ou seja, à operação, no de- 
curso da qual, sem alterar a sua estrutura, se extrai 
duma certa massa de matéria, uma fracção suposta 
representativa daquela 


1.1 — Propriedades gerais dos lotes de material a granel 
1.1.1 — Definição dum lote de material a granel 


a) Na sua definição continua o lote é conside- 
rado: 


* geometricamente: como o conjunto de pon- 
tos x interiores a um certo domínio D, do 
espaço considerado. 


* Fisicamente: como o conjunto dos valores 
tomados por um grupo de funções descri- 
tivas da matéria, presente em cada um dos 
pontos x interiores a um certo domínio D, 
do espaço considerado. 


b 


Da 


Na sua definição descontínua o lote é consi- 
derado: 


* geometricamente: como um conjunto dis- 
creto de subconjuntos contínuos. Cada sub- 
conjunto pode conter quer um fragmento 
único, quer um grupo de fragmentos. 


Eu 


Fisicamente: como o conjunto dos valores 
tomados, pada cada um dos subconjuntos 
contínuos, por um grupo de parâmetros 
descritivos da matéria constitutiva. 


1.1.2 — Expressão do teor dum lote de material a granel 
1.1.2.1 — Meio continuo 


Se a é o teor do lote L em constituinte crí- 
tico A, poderemos escrever: 


a -[, datas ff u (x) dx 
' D, | ba D, 


ou seja, é a média ponderada da função indicatriz 
a (x) — teor, em constituinte A, no ponto x, sendo a 
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função ponderada (x) — massa de sólidos por uni- 
dade de espaço. 


Por unidade de espaço entende-se: 


unidade de volume: lotes com três dimensões es- 
paciais; na 

unidade de superficie: lotes com duas dimensões 
espaciais; 

unidade de comprimento: lotes com uma dimen- 
são espacial; 

unidade de tempo: lotes com uma dimensão tem- 
poral, 


1.1.2.2 — Meio descontinuo 


(i) — Se considerarmos o lote como um conjunto 
de um número N de grupos de fragmen- 
tos G (caso dos lotes de zero dimen- 
sões) teremos: 


a = a m,/3 M, 
da tem f 
em que: 


a —teor em A do Lote L; 
a. — teor em A de G. 
M. — massa de G. 


(ii) — Se considerarmos o lote como um conjunto 
de N, fragmentos teremos: 


“= 2 M, fz M, 


em que: 


a —teor em A do lote L, 
a —teor em A do fragmento F., 
M, — massa de F 


ze 
1 


1.2 — Critérios de apreciação de uma amostragem 


Tendo a amostragem por objectivo a estimação 
de um parâmetro de qualidade desconhecido dum 
lote L e sendo ela uma operação aleatória, como tal 
geradora de erros, analisemos quais as propriedades 
de que deve gozar uma «boa» amostragem. 


1.2.1 — Propriedades definidas «à priori» 


São as qualidades que devem possuir os apare- 
lhos de amostragem e que devem ser respeitadas 
aquando da sua concepção e realização. A qualidade 
fundamental é a correcção que se define como: 


* óptica contínua: quando assegura uma densi- 
dade uniforme de probabilidade de colheita 
em todo o domínio D e nula no seu exterior. 
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* óptica descontínua: quando assegura uma pro- 
babilidade uniforme de colheita a todos os 
fragmentos que constituem o lote L. 


A amostragem diz-se incorrecta quando não asse- 
gura uma probabilidade uniforme de colheita, 


1.2.2 — Propriedades definidas à «posteriori» 


São aquelas que contam para o utilizador. 


Sejam: 

a, — teor de L em constituinte A 

E — uma amostra extraída de L 

à. — teor de E em constituinte A 

e — erro de amostragem em constituinte A 


= gd, 

Como já é sabido, na prática, a amostragem 
nunca é rigorosamente justa. Assim, consoante os 
objectivos, teremos: 


(i) — amostra comercial de rotina, entre um 
mesmo comprador e um mesmo vendedor: 
a principal qualidade deverá ser a equita- 
bilidade para o que é condição suficiente, 
mas não necessária, a amostragem ser 
justa —m (e) =0— quando o teor é 
função linear do valor comercial. 


(ii) — Amostra técnica de rotina (controlo de 
qualidade). Interessa principalmente a 
fidelidade — o (e) < 0. 

(iii) — Amostra isolada com objectivos técnicos 


ou comerciais. Interessa a precisão — 


m(ej=0e st te) e si — ou, pelo menos, 
a representatividade — m(e2) = 0º (e) + 
+m? (e) < R,2. 


Do ponto de vista prático, justica e fidelidade 
só são parcialmente dissociáveis porquanto os erros 
mais graves, que surgem quando a amostragem é 
incorrecta, afectam-nas simultaneamente. Assim, O 
objectivo fundamental deverá ser a correcção que 
constitui a única garantia de representatividade. 
Como segundo objectivo teremos a fidelidade. 


Tudo o que se segue será consagrado a determi- 
nar as condições de correcção e fidelidade, numa 
palavra, de representatividade de uma amostragem. 


1.3 — Análise dos processos de amostragem. Introdução 
aos modelos de amostragem 


1.3.0 — Todos os aparelhos e métodos de amostragem 
são baseados em um dos dois processos seguin- 
tes: processo de implantação de incrementos e 
processo de partilha 
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Ambos se encontram esquematizados na figura 1. 
Os primeiros aplicam-se a lotes não manipuláveis e 
ainda a certos tipos de lotes manipuláveis (ex.: cor- 
reias transportadoras e tomas para análise química). 
Os segundos aplicam-se exclusivamente a lotes mani- 
puláveis. 


1.3.1 — O Processo de Implantação de Incrementos 


Pode-se decompor o processo de implantação de 
incrementos num certo número de operações elemen- 
tares como segue: 


(i) — Integração: é a selecção, no domínio D, 
ocupado pelo lote L e segundo uma lei 
definida, dum certo número Q de pontos 
que constituem os incrementos pontuais. 


(ii) — «Découpe»: é a aplicação, em cada um dos 
incrementos pontuais, dum certo volume 
a que chamaremos «découpe» e que deli- 
mita um incremento modelo, 


(ii) — Toma: consiste na materialização da 
extracção dos fragmentos contidos nos 
incrementos modelo. Na toma obtêm-se os 
incrementos reais. 

(iv) — Reunião: é a operação que conduz à 
amostra real através o agrupamento dos 
incrementos reais. E uma operação exac- 
ta na sua essência, 


Destas quatro operações elementares apenas a 
integração constitui uma amostragem propriamente 
dita. As operações subsequentes limitam-se a mate- 
rializar a amostra pontual. 


1.3.1.1 — Integração 


O processo de integração insere-se na óptica 
continua já apresentada. 
Seja: 


L: Lote de material (a granel ou compacto) 
a avaliar; 

Domínio ocupado por LIL, no espaço Rº 
(p=1,2,3); 

x: Um ponto RP 

Teor no ponto x; 

Massa por unidade de espaco no ponto x; 


Um ponto de D, seleccionado, segundo 
uma certa lei de probabilidade, para 
constituir um incremento pontual; 


Q: Número de incrementos pontuais; 


E: Amostra pontual formada pela reunião de 
Q incrementos pontuais; 


a: Teor de L em constituinte crítico. 


a.: Teor de E em constituinte crítico. 
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Por definição: 


= a(x) u(x) dx (x) dx 
“E J D, : [o 


ap = a(x) n(x) / > utx) 
mM q 


A análise destas duas expressões mostra que a 
operação que consiste em considerar à, como esti- 
mador de a, se assemelha muito estreitamente com a 
estimação por pontos de um integral. 


Designa-se por lei de integração, a lei segundo 
a qual são designados os Q pontos X,. As leis mais 
usuais são: 


* integração sistemática com implantação alea- 
tória; 
* integração aleatória estratificada; 


* integração aleatória. 


O modelo de amostragem que estuda o processo 
de integração é o modelo de integração, que se apoia 
no definição continua do material a granel e na ca- 
racterização da heterogeneidade geral de distribuição. 
Poderá pois aplicar-se a todos os meios materiais, 
desde os sólidos compactos aos líquidos, passando 
pelas associações gás-sólido, líquido-sólido e líquido- 
“líquido. 

Pode ainda dizer-se que o modelo de integração 
é o estudo do problema com um certo recuo, numa 
escala tal que a estrutura fragmentária da matéria 
se desvanece, para melhor revelar a estrutura em 
grande escala da distribuição espacial ou temporal. 


1.3.1.2 — «Découpe» 


A operação pontual de integração, sendo mate- 
maticamente fácil de tratar, é materialmente irreali- 
zável. Serã necessário pois materializar os pontos 
de colheita o que é habitualmente feito em duas eta- 
pas: uma geométrica («découpe») e outra material 
(toma). 


Considere-se um certo domínio D, do espaço Rº, 
definido em forma, dimensão e orientação, e pequeno 
comparado a D a que chamaremos «découpe» (ou 
volume regularizador ou de suporte) e que é apli- 
cado a cada um dos Q pontos X, 


O incremento modelo é constituído pelo conjunto 
dos fragmentos cujo baricentro estã no interior da 
adécoupe» modelo, 

A amostra modelo resulta da reunião dos incre- 
mentos modelo. 


A «découpe» não é uma amostragem. No entanto 
é uma operação aleatória que convém representar 
por um modelo. Como ela define um grupo de frag- 
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ESQUEMATIZAÇÃO DOS DOIS PROCESSOS DE AMOSTRAGEM 
(Segundo P.Gy) 


lote 
d | 
PROCESSOS DE | | avanar PROCESSOS 
IMPLANTAÇÃO DE 
DE INCREMENTOS PARTILHA 
1) Integração — | 1) Divisão — 
incrementos pontuais Fracções — modelo 


2) “"Decoupe” — 
Incrementos modelo 


2) Separação — 
Fracções reais 


3) Escolha — 


Incrementos reais 


3) Toma — 
Incrementos reais 


4) Reunião — 
Amostra real 


4) Reunião — 
Amostra real 


Fig. 1 
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mentos, poderá ser encarada, dada a sua ambiva- 
lência, segundo duas perspectivas: 


a) Perspectiva continua 


Se considerarmos D, como o domínio regu- 
larizador, quando D,. < D,, teremos: 


À = a (x) u,(x) dx f Hotx) dx 
“4 » e / D, o 


dy — 3 a (x) nx) /z HE) 
q q 


A estimação de a,, com a ajuda de a, é uma 
integração que faz intervir as funções regu- 
larizadas, e não as funções pontuais. Como 
os variogramas experimentais são obtidos a 
partir de funções regularizadas (o suporte 
nunca é pontual) pode-se afirmar que o 
o modelo de integração regularizado cobre 
simultaneamente a integração e a «découpe». 


b) Perspectiva descontínua 


A especulação sobre a probabilidade de colhei- 
ta dos fragmentos que compõem a «découpe» 
e respectiva amostragem pode ser feita com 
a ajuda do modelo probabilista, o qual se sub- 
divide em dois submodelos: 


* modelo probabilista simples, no qual cada 
fragmento F, de L intervém individual e 
independentemente com uma certa proba- 
bilidade P de ser colhido; 


* modelo probabilista composto no qual cada 
grupo (e de fragmentos, suposto indivisivel, 
intervém individual e independentemente 
com uma certa probabilidade P de ser 
colhido. 


O problema da «découpe» situa-se a nível do 
modelo probabilista composto. 


1.3.1.3 — Toma 


E a operação material que permite extrair os 
incrementos reais. O problema do erro de toma está 
inteiramente ligado ao respeito ou desrespeito da 
regra do ressalto(!), Na prática a toma é sempre 
uma operação aleatória. 

Se a regra do ressalto fosse rigorosamente 
cumprida, a probabilidade de colheita P, de um 
fragmento contido na «découpe» seria igual a 1, e 
igual a O a probabilidade de colheita de um frag- 
mento exterior à «découpes. 


A toma é pois analisada no âmbito do modelo 
probabilista simples. 


1.3.2 — Processo de Partilha 


Os seus protótipos são a amostragem através de 
um amostrador Jones ou por pazada alternada. 

Pode-se decompor um processo de Partilha numa 
sequência de operações elementares, a saber: 


(i) — Divisão: é a demarcação do domínio 
ocupado pelo lote, através a delimitação 
de um certo número de subdomínios com- 
plementares, de extensão igual ou desi- 
gual, que constituem as fracções modelo. 

(ii) — Separação: é o isolamento das diferentes 
fracções-modelo que se transformam em 
fracções reais. 

(iii) — Escolha: é a selecção, segundo uma lei 
definida, das fracções reais que consti- 
tuirão os incrementos reais. 

(iv) — Reunião: é a operação que conduz à 
amostra real, por agrupamento dos in- 
crementos reais. 


Na realidade, a divisão é assimilável à «découpe», 
a separação à toma, a escolha à integração e a reu- 
nião à reunião. Assim, tudo o que se disse no Capi- 
tulo 1.3.1 se aplica «mutatis-mutandis» ao processo 
de Partilha. 


1.3.3 — Conclusões 


As operações elementares dos dois processos de 
amostragem estão incluidas no âmbito de dois mo- 
delos de amostragem: 


* O modelo de integração que constitui o pro- 
longamento da definição contínua da matéria, 
a granel ou não; 

* O modelo probabilista que constitui o prolon- 
gamento da definição descontínua da matéria 
a granel e que se subdivide em dois submode- 
los: o simples e o composto: 


(ii— No âmbito do modelo de integração, o 
lote L a avaliar é considerado como o 
conjunto de pontos x interiores a um 


certo domínio fechado D, C Rº. Cada 
ponto x é suficientemente caracterizado 
pelos valores de duas funções: a(x) teor 
e u(x) massa por unidade de espaço (ou 
de tempo); 

(ii) — No âmbito do modelo probabilista sim- 
ples, o lote L, é considerado como um 
conjunto discreto de N/ fragmentos F.. 
Cada fragmento F, de L é suficiente- 


(1) Regra segundo a qual os fragmentos ao chocarem sobre uma aresta horizontal ressaltam para o lado em que 


ge encontra o seu baricentro. 
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mente caracterizado pelo valor de dols 
parâmetros: a sua massa M, e o seu 


teor d,) 


(iii) — No âmbito do modelo probabilista com- 
posto, o lote L, é considerado como um 
conjunto discreto de N grupos de frag- 
mentos G|, supostos indivisíveis. Cada 
grupo G é suficientemente caracterizado 
pelos valores de dois parâmetros: sua 
massa M eo seu teor a. 


2. HETEROGENEIDADE DE UM LOTE DE MATERIAL 
A GRANEL 


2.0 — Introdução 


A amostragem dum material rigorosamente 
homogéneo é uma operação exacta. Os erros de 
amostragem são devidos à heterogeneidade da ma- 
téria e a sua multiplicidade está logicamente rela- 
cionada com a multiplicidade de formas de que a 
heterogeneidade se pode revestir. 

Daí ser lógico começar por aprofundar a noção 
de heterogeneidade antes de desenvolvermos os mo- 
delos de amostragem. 


Sr 


| 
2.14 — «Aproach» qualitativo da noção de heterogeneidade 


* Diz-se que um lote de material contínuo (a 
granel) goza de constituição homogénea quan- 
do todos os seus elementos constitutivos (frag- 
mentos que o compõem) são idênticos entre 
si. Caso contrário, diz-se que possui uma cons- 
tituição heterogénea. 


* Diz-se que um lote de material, contínuo ou a 
granel, goza de distribuição homogénea quan- 
do um volume corrente, de dimensões cons- 
tantes, levado a ocupar todas as posições 
possíveis no espaço do lote, possui uma com- 
posição média constante. 


Pode pois dizer-se: 


— A heterogeneidade de constituição é uma ca- 
racterística intrínseca e inalterável da matéria, sobre 
a qual a homogeneização não tem qualquer influência, 


— A homogeneidade de distribuição é função da 
heterogeneidade de constituição, da distribuição espa- 
cial dos fragmentos e da escala de observação, ou 
seja, da dimensão do volume corrente. O seu limite 
superior é a heterogeneidade de constituição, 


(!) Para mais detalhes vide Bibl. 22 (pgs. 16 a 27). 
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2.2 — Classificação e propriedades das distribuições de 
fragmentos 


2.2.4 — Introdução 


Embora nos interesse mais particularmente a 
função a(x), as conclusões são também válidas para 
a função q(x). 

O espaço no qual existem os lotes está mergu- 
lhado no campo gravítico, pelo que é inevitável que 
a segregação gravítica tenda a classificar vertical- 
mente os fragmentos de acordo com a sua dimensão, 
massa volúmica, etc. sendo a principal responsável 
pela alteração da homogeneidade de distribuição. 

A classificação que se indica parte da ordem 
perfeita, que corresponde a uma segregação com- 
pleta e a uma relação biunívoca entre a personali- 
dade de um fragmento e a sua posição no lote e 
chega à desordem perfeita que corresponde a uma 
correlação nula entre a personalidade de um frag- 
mento e a sua posição. 

Os casos práticos estão obviamente situados 
entre estes dois limites e são ('): 

* Desordem local sobreposta a uma ordem pro- 

funda; 

* Ordem relativa; 

Desordem local sobreposta a uma ordem rela- 
tiva (distribuição ordenada); 

* Ordem relativa delineada e redistribuiída; 
Desordem local sobreposta a uma ordem rela- 
tiva delineada e redistribuída (distribuição de- 
sordenada ); 

* Distribuições periódicas. 


2.2.2 — Expressão mais geral de a(x) 


No caso mais geral pode-se decompor a fun- 
ção a(x) numa soma de quatro termos: 


a(x) =a, +a,(x) +a (x) + a,(x) 


em que: 


a | — constante e igual ao teor médio do lote L, 

a, (x) — componente não periódico de ordem (ter- 
mo regionalizado) 

a (x) — componente de desordem (termo local) 

a.(x) — componente periódico de ordem (termo 
periódico) 


Alguns destes componentes podem ser nulos. 
No quadro 1 indicam-se resumidamente as carac- 
terísticas das principais distribuições. 


2.2.3 — Hipótese ergódica e estacionária 


As hipóteses ergódica (similitude entre a média 
temporal e a média estatística) e de estacionaridade 
da média e da variância só poderão ser consideradas 
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válidas para a(x), quando o lote tiver sido submetido 
a uma tentativa de homogeneização. Contudo, são 
perfeitamente aceitáveis para a função a(t) quando 
se trata de um tal-qual de mina, um concentrado ou 
estéril duma lavaria, uma alimentação duma meta- 
lurgia, um minério comercial, da alimentação duma. 
fábrica de cimento, de vidro ou de cerâmica, etc. 


2.3 — Caracterização da heterogeneidade geral de dis- 
tribuição 


2.3.1 — Definição do variograma integral 


G. Matheron mostrou que era possível caracte- 
rizar uma distribuição através da função intrínseca 


QUADRO 1 


RESUMO DAS DISTRIBUIÇÕES 


ORDENADA 
j Completamente Completamente 
Distribuição avitenaão Desordenada desordenada 
Não periódica Periódica 
Não nulo Não nulo Nulo ou não nulo 
a (x) Identicamente nulo 
A monótono não monótono e não monótono 
a (x) Qualquer Mínimo 
a,(x) Identicamente nulo Não nulo Identicamente nulo 


2.2.4 — Propriedades dos componentes de a(t) 


a, (t) — componente estocástico a longo termo. 
Traduz as tendências profundas da a(t) 
bem como a continuidade das variações 
a grande escala, A sua curva represen- 
tativa é contínua e desenvolve-se em 
torno de constante a. É integrável e 
independente do suporte Do. 


a tt) — componente estocástico a curto termo. 
Traduz os fenómenos locais, as perturba- 
ções a pequena escala e a natureza dis- 
creta da matéria. Representa-se por uma 
linha quebrada em torno de a(t) +a (t). 
É integrável e depende do suporte De. 


a,(t) — componente periódica. Pode pôr-se sob a 
forma duma soma de teores do tipo 


a(t)=a, sen wt+a”. cos mwt 


E contínua, integrável e independente 
de D,. 


à 


ou semi-variograma e a partir dele deduzir a disper- 
são dos erros de amostragem. 


Seja: 

L — Lote de minério a granel em desloca- 
mento 

D, — A duração do deslocamento de L (*) 

8 — Intervalo de tempo 

a(t) — O teor do material que se desloca no 


instante t. 

3(tg)— O acréscimo algébrico da função a(t) 
entre os instantes t— 9/2 e t+ 9/2 
(to) =a(t+ 9/2) — a(t— 8/2) 

V(9) — O semi-variograma integral de a(t) 


- 


vid da 8º (t,6) dt 
D, —»+00 2D, JD 


E 


(tt) Supor-se-á a origem dos tempos coincidente com a origem do deslocamento. 
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2.3.2 — Propriedades do variograma integral (') 
(i) — Atendendo a que 
a(t) — a + a(t) + alt) + alt) 
teremos: 


VO=S VD +HS D Vosr(o) 
i : w 
Se considerarmos que a(t) é estacionária 
no 2.º momento e ergódica e que, em cada 
instante, as realizações dos componentes 
5, (t,6) são independentes, podemos con- 
cluir que V,, (6) =0e que 


V(g)=5 Vo) 


(1) — Se for: p, (8) O coeficiente de autocorrela- 
ção à escala do intervalo g e q? a espe- 
rança matemática de a 2 tj), teremos: 
Para j=1 e 2, teremos; 

lim  p;(9)=0 
a <l E— 00 


O < Vo) <o;? lim V;(0)=02 


A+ 09 


(iii) — Propriedades de V'(9) e seus componentes. 
O variograma V (4) e seus componentes 
gozam das seguintes propriedades: 


1)—0< Vi(6) <a; lim Vi(0)=0 ; 


6—00 V.lg)=g; 
tim. 
2)— MW 0 De ; Vilo)=0; 
3) — V;lg)—0;A1-COS (y 6) 
4) — V(g)=V,(6) + Vito) + Vlo) 


lim V(g)=02; Vo> 0,=V(g)=a; + 
O—E 


+o2+o; (1-cos 6) 


Nas figs. 2 e 3 estão indicados alguns exemplos 
de variogramas (segundo P. Gy). 


2.3.3 — Determinação experimental do variograma 
Distinguiremos as seguintes formas do varlo- 
grama: 
(i) — Variograma integral V(g). É a função 
definida acima; 
(ii) — Variograma discreto V'(j9,). E o estima- 
dor pontual do variograma integral(2); 


VARIOGRAMA V(9) E SEUS COMPONENTES 
CASO GERAL 


Er a “portêe 


uol [3]: w [8] v [e] 


=." 


T = periodo 


[2] : v; [8] 


[4] : v [e] 


v[e] = wlo] + w[0] + v [0] 


2 
[a] = 0, + 0,+0; 


Fig. & 


(y Vide Bibl. 11, pgs. 124 a 163 e Bibl, 22, pg. 28. 


(2) De aplicação imediata nos lotes com zero dimensões. 
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VARIOGRAMA Ví 


vIe) (1): VIO] = Voo] 


(Efeito de 


Pepita Puro) 


)- casos particulares 


vI691] = v,/ 0)» 


[2] 


V316] 


E + 
, 
N / 
N / 
EE = Ne! 
2 
+ o, 
o a 
[3] V([60] =V, [0]+ v,i6] E 
he o pa a o vie) [6] - V(0) = V,(8) + Vote) 
E ão al iIvarioagrama nado neriôdico ) 
f 
/ / Gp. === =—— RCas a = 
/ 
/ / 
/ 
a / 
) 
/ 
À / 
A 
2 
e) 
| o? = oi « oi 
2 2 . o? 
0“ mo; + 0; 
" e) 
E o E a 
Fig. 3 


(iii) — Variograma experimental v” (Jo,)). E o 
estimador experimental bruto do varlo- 
grama discreto; 

(iv) — Variograma corrigido V (J9,). E o vario- 
grama experimental após correcção dos 
erros de determinação; 

(v) — Variograma interpolado V (9). E a rela- 
ção funcional que é possível deduzir dos 
pontos conhecidos do variograma corri- 
gido V (j9,) e que constitui a represen- 
tação do variograma integral, a partir da 
qual poderemos calcular os momentos do 
erro de amostragem de a(t). 


A cu 
2.3.3.1 — Estimação de 0;2=V;(8) 


1. Cálculo de V, (6)=0,2=% (6,) 


2. O patamar do variograma ga. 


V)696,>06,: VIoy)—Vw)>O 


(1) No caso em que V, (8) %0, remete-se o leitor para a bibliografia 


2.3.3.2 — Estimação do variograma interpolado 


Apenas se indica o caso em que o termo perlió- 
dico é nulo (1). 


(i) — Cálculo de V (jo) =V(jo)—V(6) 

(ii) — A partir dos pontos Vo j9,) assim defi- 
nidos e de V (9) = 0, deduz-se a compo- 
nente V (6) do variograma interpolado. 


E habitualmente possível adoptar uma 
das expressões seguintes: 


Vl) =V,o 
— A 
Vi()=V, 
Vil) =V,0+V,62 


cujos parâmetros podem ser determina- 
dos quer por método gráfico quer pelo 
método dos mínimos quadrados. 


(11) — Como V, ()=V (9,), teremos: 


Vi) =V,(6) +V, (6) 


(Bib. 11, pes. 140-145). Também em Bib. 11, 


pes. 145-161, se descreve a organização de uma experiência variográfica e um exemplo de análise completa dum variograma 


periódico. 
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2.4 — Caracterização da heterogeneidade local de distri- 
buição 
2.4.1 — Variograma das distribuições não ordenadas 
São caracterizadas por a (x) = a (x) =0, 
donde a(x) = a, +a, (x) 


sm 


| [a,(x)]º dx= constante 
D,, 


1 
e Vi)=V lo)=0,Ã= 5 
L. 


sá DO rum 


2.4.2 — Análise da noção de heterogeneidade duma dis- 
tribuição não ordenada (") 


24.21 — Introdução e expressão matemática das dife- 
rentes formas de heterogeneidade 


Considere-se um lote de material a granel (ex.: 
minério). Para comparar as diferentes formas de 
heterogeneidade teremos que nos fixar sobre a popu- 
lação dos fragmentos isolados e sobre as populações 
dos grupos de fragmentos. Seguidamente utiliza-se 
a técnica da análise da variância e caracteres adi- 
tivos (volumes ou desiquilíbrios relativos das massas 
do constituinte crítico). 

Chega-se assim à expressão matemática das 
diferentes formas de heterogeneidade: 


* heterogeneidade de constituição do Lote L;: 


1 
Et) = (d,) me 3 3 (49 — 4 
n Í 


É uma característica intrínseca do estado de 
fragementação do lote; 
* heterogeneidade de constituição do Grupo G, 


1 


Ho (G) = 2 AA 


n 


* heterogeneidade média de constituição dos 


grupos G.. 
E a média de H, (G,) ponderada pelo seu 


efectivo N. 


Es 1 1 
Rot) = 2" AG) = pá 2 (A, —A,)? 
a J 


* heterogeneidade de distribuição de L entre os 


os grupos G, 
O > N 4 
“LL n 
em que: 


da UA == ã,) M, / 4, M,. 
An; = (Mn; — a) My, /8 M,. 
N| — efectivo do grupo G.. 


n 
N, — efectivo do lote L. 


M 


L 


& teor do lote. 


a) —teor do fragmento j do grupo G.. 


— massa do fragmento médio do lote. 


a —teor do grupo G. 
M, — massa do grupo G,. 


n 


M, , — massa do fragmento j do grupo G,. 


com a seguinte relação entre elas: 
HAL) =H(G) +H(G,) 


2.4.2.2 — Homogeneidade natural de distribuição. Parâme- 
tro de agrupamento 


Diz-se que um lote apresenta uma distribuição 
naturalmente homógena quando há independência 
entre a personalidade de um fragmento e a posição 
do seu baricentro. A distribuição espacial dos frag- 
mentos é uniforme ou, por outras palavras, os frag- 
mentos estão distribuídos aleatoriamente pelos grupos 
G, que se podem formar no seio do lote L, 

Consequências: pode-se considerar o lote L como 
um dos lotes possíveis de efectivo N, extraído ao 
acaso e um a um duma população Lo de efectivo 
infinito. 

Seja 4, a heterogeneidade do fragmento F, de L, 
e o? a variância de à, nesta população de efectivo 
infinito. 

A hipótese da homogeneidade natural é equiva- 
lente à colheita ao acaso e um a um dos fragmentos 
constituintes dos grupos G,. 

Sabe-se da estatística clássica que: 


(D— 3 3 (AA [2 =N Ho (L) 70º 
no 


se distribui em x? com N,—1 graus de 
liberdade 


D— 3 > (4402/02 =N Ho (Go) /02 
E] J 


se distribui em x? com N,—N graus de 
liberdade 


(iii) — 2 Na (4,—8,)2/0º =N,H, (G)/0 
n 


se distribui em x? com N-1 graus de 
liberdade. 


Por outras palavras: 
(i)—a variância total N//N, |, [Ho (L)], 

(ii) —a variância residual N,/N, , [H.(G)], 
(iii) — a variância entre estratos N, /. ,[H (G)| 


(1) Dado este assunto já ter sido abundantemente tratado, limitar-nos-emos aos passos e conclusões essenciais, reme- 


tendo o leitor interessado para a bibliografia. 
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constituem, três estimadores independentes e não en- 
viezados da mesma variância g2, 

Como H, (L) é um dado do problema, pode-se 
dizer que a homogeneidade natural de distribuição 
é caracterizada por: 


n[HL(G) "2 m ssmame=* ma 


ud Ed 


N-1 1 | 
ni ]= a Pott)=Ho(G, mino + Hot 
L l 


O parâmetro de agrupamento y, característico 
da escala de «découpe» (do suporte) será definido por 


v = (N—N) /(N-d) 
2.4.2.3 — Caracterização do grau de segregação 


No domínio natural, pode-se escrever; 
| Ii 
HaL) > KH (G,) a Lry "o (L.) 


O parâmetro de segregação £ característico do 
grau de heterogeneidade da distribuição é definido 
por: 


EI = Up A) 
£ pode tomar os seguintes valores: 


£=1—a distribuição é rigorosamente hete- 
rogênea 


O < E< 1-—heterogeneidade natural da distri- 
buição 


é — 0 — homogeneidade natural 
—1/y < & < 0 — hiper-homogeneidade (artificial) 


É = — 1/y — homogeneidade rigorosa 


Esta expressão contém efectivamente os três 
elementos responsáveis pela heterogeneidade de dis- 
tribuição: 


H, (L) — heterogeneidade de constituição 


É — segregação 


dimensão dos grupos, ou escala de 
observação. 


F 


Os domínios de H (G,), m(G,) e é estão esque- 
matizados no quadro 2. 

A influência da escala de observação está esque- 
matizada no quadro 3. 


2.4.2.4 — Expressão geral de Hp (Gn) 


Explicitando a expressão do capítulo anterior 
em função das características do fragmento genérico 
F. de L, obtém-se a seguinte expressão geral da he- 
terogeneidade de distribuição do lote L, que será utl- 
lizada na teoria de amostragem: 


Er No, " a-a, M, . 
I(G)=(1+y O) (N-1) — A ( es 
À a 


N,—1 a, M, 


QUADRO 2 


DOMINIOS DE H, (Gn), Ho (Gn) e £ 


Heterogeneidade Heterogeneidade 
natural de 


distribuição 


rigorosa de 
distribuição 


d-f=y 


Homogeneidade 
natural de 
distribuição 


Hiper-homogeneidade Homogeneidade 
(artificial) de 


distribuição 


rigorosa de 
distribuição 


H, (L) 


Distribuição Distrito Distribuição | 
perfeitamente na ç perfeitamente Distribuições artificiais 
ordenada desordenada 1 Rag fr sem existência natural 
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QUADRO 3 


INFLUÊNCIA DA ESCALA DE OBSERVAÇÃO. 
VALORES LIMITES DE Y Hp (Gn) e Ho (Gn) 
NO DOMÍNIO O < £ <1 


| ME E a ESA a * 
N, > N > 1 
à < Y < to 
H, (L) o Ei) s 0 
0 < H(G) e H(L) 


3. PROCESSOS DE IMPLANTAÇÃO DE INCREMENTOS 
— MODELO DE INTEGRAÇÃO 


3.0 — Introdução 


O modelo de integração resolve o problema da 
amostragem na perspectiva contínua, o que lhe dá 
uma grande generalidade. 

Os resultados do estudo deste modelo aplicam-se 
indiferentemente a todos os tipos de teores e nomea- 
damente: 


* teores em componentes químicos ou mine- 
ralógicos. 

* teores de humidade, 

* proporções das classes granulométricas. 


Na sua essência o problema é rigorosamente o 
mesmo, qualquer que seja o número de dimensões 
do espaço do lote. Para fixar ideias estudou-se o caso 
importante do espaço com uma dimensão temporal 
que serve de modelo aos lotes em deslocamento (1). 


3.1 — Leis usuais de integração 


Às leis mais usualmente utilizadas para designar, 
no seio do domínio D, ocupado pelo lote L, os pontos 
sede dos incrementos pontuais são, como se sabe: 


* integração sistemática com implantação 
aleatória. 

* integração aleatória estratificada. 

* integração aleatória. 


Podemos, resumidamente, dizer o que se passa 
quanto às três leis: 


(i) — Integração sistemática e aleatória estati- 
ficada. 


a) Quando. D /9 é um número inteiro 
qualquer, a densidade de probabilidade 
de implantação g(x) é uniforme e igual 
a i /8 em todo o domínio D,. O núme- 
ro Q de incrementos é constante e 
igual a D,/0. 


b) Quando D, /9'é um número não intei- 
ro > 50, a densidade de probabilidade 
de implantação g(x) é praticamente 
uniforme, embora não o seja de modo 
rigoroso, O número Q de incrementos 
é uma variável aleatória podendo to- 
mar os valores Q' e Q+1. 


c) Quando D,/s é um pequeno número 
não inteiro, a densidade de probabili- 
dade de implantação g(x) já não é 
sequer praticamente uniforme. O nú- 
mero Q de incrementos é uma variável 
aleatória podendo tomar os valores Q' 
e Q'+1. 


dj Só são correctas se e só se g for um 
submúltiplo de D,. 


(ii) — Integração aleatória. 


O número Q de incrementos é o 
parâmetro determinante da lei e é cons- 
tante. Donde, a densidade de probabili- 
dade de implantação g(x) é uniforme 
em todo o domínio D, e igual a Q/D,. 
A integração aleatória é sempre correc- 
ta. 


3.2 — Desenvolvimento do modelo de integração 


3.2.1 — Distribuição da massa da Amostra Mp 


(i) — No caso que nos ocupa, o variograma 
Vy(9) pode, na prática, ser reduzido, 
numa região limitada, a uma expressão 
linear da forma 


V(6) = Vw 8 + Vw 


(ii) — Na hipótese de Q > 50, a massa M, da 
amostra. E segue uma lei praticamente 
normal, cujos momentos são: 


(1) A amostragem primária dos lotes a 1, 2 ou 3 dimensões espaciais deve ser resolvida com base nos trabalhos de 


Matheron. 
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* Integração sistemática: 


Vu Vuz 
m (Mp) =D, Dc/9; (Mp) =D |—— + — 


* Integração aleatória estratificada: 


Vu Vus 
m (Mp) = M, D./0; o AM,) = D, 3 +- 


* Integração aleatória: 


j Vm D, 
mM) =M D./D,; (Mp) = (—— + Vi) 


equivalente à da estatística clássica 


2UM,) = Q cXM,) 


3.2.2 — Distribuição de Ap — Correlação Ap, My 


(ij) — O que se disse para a distribuição da 
massa da amostra é válido para a distri- 
buição da massa do constituinte crítico A. 

(ii) — Quando as componentes periódicas a,(X) 
e ua(x) são nulas, é habitualmente possível 
admitir como hipótese de trabalho que os 
incrementos G e G , duma mesma amos- 
tra E têm características independentes. 
Pode pois escrever-se 


p(Ã Me) =e(Ap MC) 


3.2.3 — Distribuição de ar 


3.2.3.1 — Integração correcta. Caso geral 
Def.: Wx E D, ; go) =E, 
IJ xe D: uíx) Eu 
p(uíx), a(x)) £o 
a(x) = a, 


(1) — Qualquer que seja a lei de integração, a 
variância relativa e o teor de a, podem 
exprimir-se, em 1.º aproximação, a partir 
das variâncias relativas de M e A, e da 
correlação p(A,M) calculáveis a partir 
duma experiência variográfica 


UZ(a,), = UZa,) + UMM,) —2 p(A,M). U(AÇ). U(M,) 


Esta primeira aproximação é sempre su- 
ficiente na prática(2). 


(iij — Quando a lei de integração é correcta, a 
integração é justa. 


(ili) — Qualquer que seja a lei de integração, 
desde que assegure densidade uniforme de 
probabilidade, a integração é justa em 2.º 
aproximação unicamente nos 3 casos se- 


guintes: 
u(x) = cte 
p(alx), u(x)) — 0 
a(x) = cte 


Em todos os outros casos é enviezada, podendo 
o viés ser calculado a partir de uma experiência 
variográfica. 

A origem profunda deste viés reside na natureza 
da função teor. Um teor médio é uma média ponde- 
rada e não aritmética, excepto nos 3 casos apresen- 
tados, 


3.2.3.2 — Integração correcta. Ponderação constante 


W XE D,: P(X) — E, 
p(x) = u, => Vy(9) = (distribuição homogénea) 


* média: 


Quando a função de ponderação (x) e a densi- 
dade de probabilidade g(x) são uniformes no domi- 
nio D,, o teor a, segue uma lei normal de média igual 
à Bs 
A integração é rigorosamente justa, 


* Variância. 


As expressões da variância apresentadas no caso 
geral continuam válidas(3). Basta notar que: 


6) = 2, V.(0) 


A 


< 


3.2.3.3 — Integração correcta. Correlação nula 


Wal gix)=8 
p(a(x), uix)) — o 
O teor a, segue uma lei praticamente normal de mé- 


dia m(a,.) =a e de variância relativa 


UZ(a,) = UMA,) — UMM,) 


3.2.3.4 — Integração qualquer. Teor constante 


v(x) e D: a(x) =a 
ap=a, m(a,) =a, 
o(a) = o 


(2) Para explicitação em função das três leis de integração e para o cálculo da 2.º aproximação, vide bibliografia 


—Bit. 11, pãg. 213-217. 
(3) Vide Bibliografia Bib. 11 —pág, 219. 
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3.2.3.5 — Integração incorrecta. Ponderação constante(!) 


YxeD ; gix) E, 
(x) = 


Na integração sistemática e aleatória estratifi- 
cada quando 6 -—> 0 e Q' é grande (> 50) pode admi- 
tir-se que o desvio sistemático é desprezável, isto é, 
a integração é praticamente justa. 

A integração aleatória é sempre justa, 


A bibliografia (Bib. 11, pág. 229-241) apresenta 
numerosos exemplos de aplicação. 


3.3 — Critérios de apreciação de uma amostragem na 
óptica do modelo de integração 


Como já vimos a hierarquia das qualidades glo- 
bais de uma amostragem é: 


1.º exactidão, 
2.º precisão. 
3." representatividade. 


As condições práticas são: 
(i) — Integração exacta. 
Para que a integração seja exacta é 
necessário e suficiente que a matéria a 


amostrar tenha constituição rigorosamen- 
te homogénea: 


VOU) eD =a(x) = ay 


No caso geral, a integração nunca é 
exacta, 


(il) — Integração precisa. 


É necessário: 


a) que seja rigorosamente justa, o que 
implica que a densidade de probabili- 
dade de implantacão de incrementos 
seja uniforme: 

WxeD —>gix) =£&. 
WYxED —>g(x)=o0 


b) que seja satisfeita uma das duas co0n- 
dições suficientes seguintes: 


b 1) WYxeD >ulx) =y 
b. 1) YxeD, — pla(x), u(x)] = o 


c) que seja fiel, isto é, que a lei de inte- 
gração e o seu parâmetro determinante 
sejam escolhidos em função do nível 
de fidelidade fixado. 


(1) Para mais detalhe vide Bib. 11, pág. 223-229, 


(iii) — Integração representativa. 
É necessário: 


a) que seja justa o que implica a correc- 
cão(2) 


YxeD >g(x) =g, 
Yx$éD —g(x)=o 


b) que seja fiel (vide (il) — c). 


34 — Análise do erro total de integração 


Tudo se passa como se o erro total de integração 
EIl= (a —a, )/a fosse a soma de dois erros inde- 
pendentes: 

* O erro de integração da função 

(x) = EI, = (ap—a, )/a, 
em que a, é uma média aritmética. 


* O erro da função ponderadora u(x) — ED. 
ED é definido de tal modo que 


m(EI) = m(ED) + m(EL) 
(EI) = c2(ED) + cHEI ) (3) 


Por outro lado, o erro da função a (x) — SE. = 
pode ser considerado como a soma de três erros: 

* Erro de integração do termo regionalizado a (x) 
local a. (x) 
periódico a,(x) 


+ > » » >» % 


. >» » > + > 


Us dois primeiros momentos destes erros, supos- 
tos independentes por hipótese, são: 


m(EI ) = m(EL) + m(EL,) + m(EL,) 
(EI) = cXEI,) + c2(EI,) + cXEL,) 


O erro total de integração pode pois considerar- 


-se como a soma de quatro erros independentes, o 
que permite somar os seus dois primeiros momentos. 


3.5 — Erro de Ponderação ED 
(i) — Média. 
m(ED) = m(EI) 
Isto é, o viés do erro total de integração é 


inteiramente devido às flutuações de (x), 
só se anulando quando (x) = constante. 


(2) Esta condição é satisfeita para as três leis de integração, desde que D,/6 seja inteiro ou maior que 50. 


(3) HEI) = Ulag). 
oHEI) = Uaç'). 
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Daí o interesse em regularizar o cau- 
dal de forma precisa, antes da amostra- 


gem. 
(ii) — Variância, 
Por definição c(ED) = UZa,) — UZa,'). 


Quer a média, quer a variância de ED podem 
ser determinadas através da realização duma expe- 
riência variográfica. Para os leitores interessados 
recomenda-se Bib. 11, pág. 265-275. 

Na prática, os erros de caudal cometidos, quando 
este é regularizado, ponderal ou volumicamente de 
forma precisa, são aceitáveis. 


3.6 — Erro de integração do termo regionalizado a;(x) 


4.6.1 — Definição 


Chama-se termo regionalizado à componente 
a,(x) da função a(x) que tem por variograma a com- 
ponente V,,(0) de V (8). Pode-se determinar a com- 
ponente V,,(6) do variograma interpolado a partir 
do resultado duma experiência variográfica. 

Ao longo deste capítulo, consideram-se as se- 
guintes hipóteses de base: 


(i) — x pertence a um espaço com uma dimen- 
são temporal; 

(11) — uíx) = cte = Hr 

(ii) —a(x)=a +a (x) ; Va(0)=V(0); 


(iv) — A densidade de probabilidade de implan- 
tação dos incrementos é uniforme. 


YxeD >gix) = 8 
x é D —glx) o 


(v) — A componente regionalizada do variogra- 
ma é linear, no domínio que nos interessa, 
admitindo por equação 

VO<D : V, (6) =V,º 
expressão, na qual, V,, é o parâmetro de 


regionalização linear. Esta hipótese é 
confirmada, com frequência, na prática. 


3.6.2 — Média do erro El, 
Devido a (ii) e (iv): 


m(EI,) = 0 ou m(a,) = 4,1 


it) Vide bibliografia, 
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3.6.3 — Variância do erro El; 
Como a (x)=o e V.=0, teremos: 


(1) — Integração sistemática 


GABI) — —e— = — 
O QigLa 6açiDy 
(li) — Integração aleatória estratificada 
Va, D. Va, b* 
(lii) — Integração aleatória 
Va, Dr, 


(EL) =——— 
30a,* 


Uma experiência variográfica permite determi- 
nar Va 

Verifica-se pois que a integração sistemática é 
duas vezes mais fiel que a aleatória estratificada e 
2 Q vezes mais fiel que a aleatória. 


* Condições de anulação de (EI). 
Só se anula quando Va, o que corresponde 
ao caso particular a, (x) = &,, — constante. 


* Condições de minimização de c(EI,). 


q) Escolhendo a lei de integração mais fiel com- 
patível com os outros dados do problema, 

b) Escolhendo o mais pequeno valor de 8 com- 
patível com os outros dados do problema. 


3.6.4 — Exemplo de aplicação 


Tomemos o exemplo de uma experiência vario- 
gráfica efectuada por P. Gy(') e consideremos uni- 
camente os dados pertinentes, 


D, = 480 min (um turno de 8h): 


a = 1,58 X 10% (g Pb/g sólidos); 
v6/o<60<140min: V,(6) = 


= 0.078 x 10º 6 
a) Cálculo das variâncias, 


Das expressões indicadas em 3.6.3 calculam-se 
imediatamente as variâncias: 


o? = 1.085 X 10/62 
02 = 2.17 x 107/62 
q; — 0.05/0 
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b) Cálculo do intervalo 8 do comprimento do es- 
trato ou do número de incrementos necessários para 
que a variância relativa seja < q). 

De 3.6.3 — tem-se imediatamente: 

(1) — Integração sistemática. 


9<0,=0 à (6 D//V,)*=3040 q, 


(ii) — Integração aleatória estratificada. 


<= & (3 Di /Va)” = 2150 6, 


(iii) — Integração aleatória. 
Q 2 o e ia D,/3 Co? a? = 0,05 Co 


Assim, supondo que se desejava fixar um inter- 
valo de confiança de + 1% (ao nível 25) ter-se-á 


que tomar para = 010”, o que permite calcular: 


i) 15.2min ou 95 32; 
ii) 10.75 min ou Q > 45; 
ii) Q > 2000. 


3.7 — Erro de integração do termo local a(X) 
3.71 — Definição 


Chama-se termo local à componente a, (x) de 
a(lx) que admite por variograma a componente 
Va(0) de V,(0). Pode-se determinar a componente 
Vas(B) = V,, do variograma interpolado à custa dos 
resultados de uma experiência varlográfica. 

Ao longo deste capítulo fazem-se as seguintes 
hipóteses: 

ij) x pertence a um espaço com uma dimensão 
temporal; 

ii) u(x)=p =cte; 

HD a(x) =a, +ax); 


donde: VA = Vo t) = Voy expressão na qual V,, 
é o parâmetro de desordem local. 


3.7.2 — Média do erro El; 

Como (x) =p, e como as três leis da integração 
asseguram uma densidade umiforme de implantação 
dos incrementos 


m(EI,))=o ou m(as )=a,, 


desde que D /0 seja inteiro ou maior que 50. 


3.7.3 — Variância do erro El; 
Pela aplicação das leis de 3.2.3.2, tem-se 
z e e 2” 2 
cj(EL,) =2;(E1,) =c,(EL,) EM ui 9/D, a =V.,/Q E 


Uma experiência varlográfica permite determi- 
nar Va. 


* Condições de anulação. 
Só se anula quando V,=0(a,(x)=a,,)(1) 


* Condições de minimização. 


Escolhendo o menor valor de 6 compatível com 
os outros dados do problema. 


3.7.4 — Exemplo de aplicação 
Retomemos o exemplo referido em 3.6.4. 


D, =480 min 
a =1.58xX10º (g P,/£ sólidos) 
= O 


(2) 
Ve=Vo-o? =10-0.25x10 *=0.75x10 


(i) — Cálculo da variância. 


De 3.7.5 tem-se 


q (EL) =6.25x 10º 8 = 3000xX10* /Q 


(ii) — Se se fixar uma tolerância relativa 
+20,=*H2 % (5? =10"*), ter se-á: 
0 <0,=o; à D /V,=16X10º 9;=16 min 
Q>0,=V.,/028/=3000X 10º /0,=30 


Verifica-se facilmente que Q,X8,=D,. 


3.8 — Erro de integração do termo periódico a;(X) 


Em «Théorie et Pratique de INEchantilonnage 
des Matiêres Morcelées» de P. Gy, trata-se desenvol- 
vidamente o tema na hipótese da amplitude e fre- 
quência de a, (X) serem constantes, 

Como mem sempre isso sucede e como é por ve- 
zes muito difícil calcular o período do fenómeno sem 
recorrer a experiências variográficas longas e onero- 
sas, apresentam-se as conclusões e recomendações 
tendentes a diminuir o valor da variância do erro da 
integração do termo periódico. 


() Veremos no âmbito do modelo probabilista, quando se desce à escala do fragmento, as condições em que esta 


variância é susceptível de se anular. 
() Variância dos erros de preparação e análise, 
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(i) — Quando é impossível definir com precisão 
o período T do fenômeno periódico e quan- 
do o risco representado pelo valor máximo 
da variância 


% (0/22 


é inaceitável, deve evitar-se a integração 
gistemática. Escolher-se-á uma amostra- 
gem aleatória estratificada cujo valor máã- 
ximo da variância é Q vezes menor. 

A integração aleatória nunca é me- 
lhor que a aleatória estratificada. 


(ii) — Quando é possível definir com precisão o 
período T, pode usar-se a integração sis- 
temática desde que 6 seja um submúltiplo 
de T ou um múltiplo ímpar de T/2. Tam- 
bém a amostragem aleatória estratificada 
conduz, em média, a uma variância me- 
nor. A integração aleatória nunca é me- 
lhor que a aleatória estratificada. 


(iii) — Quando 14 (a,/a )2 é aceitável, pode 
usar-se a integração sistemática sem pe- 
rigo. 

(iv) — Hã sempre interesse em procurar e tentar 
eliminar, se possível, as causas das varia- 
ções periódicas da qualidade da matéria 
a amostrar, devido aos riscos inerentes 
à sua amostragem. 


(v) — É conveniente não esquecer que se integra 
simultaneamente os termos a (X) e a (X). 
O que importa é minimizar a soma 
c(EI,) +oZ(EI,) utilizando, segundo os 
casos, a integração sistemática ou a alea- 
tória estratificada. A integração aleatória 
conduz sempre a maior variância total. 


(vi) — Na dúvida, ou na falta de elementos pre- 
cisos, a integração aleatória estratificada, 
embora não seja sempre a mais fiel, cons- 
titui sempre a melhor solução de compro- 
misso conduzindo a amostragens nem de- 
masiadamente imprecisas nem demasiado 
caras. 


4. PROCESSOS DE IMPLANTAÇÃO DE INCREMENTOS 
O MODELO PROBABILISTA 
4.0 — Introdução 


O modelo probabilista resolve o problema da 
amostragem na perspectiva descontínua. 

Os resultados aplicam-se indiferentemente a to- 
dos os tipos de teores e nomeadamente: 


* teores em componentes químicos ou mine- 
ralógicos; 

* teores de humidade; 

* proporções das classes granulométricas, 
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Resolve os problemas nos quais intervêm quer os 
fragmentos isolados, quer os pequenos grupos de 
fragmentos que constituem os incrementos. São eles: 
a análise do erro de integração do termo local a, (x), 
a «découpe» e a toma. 


4.1 — Desenvolvimento do modelo probabilista 
41.1 — Definição 
Distinguiremos dois sub-modelos: 


(1) — Modelo Probabilista Simples: considera o 
lote L, como o conjunto dos N, fragmentos 
F. que o compõem. Cada fragmento parti- 
cipa indvidual e independentemente na 
amostragem com probabilidade P, bem 
definida, de ser colhido. 

(ii) — Modelo Probabilista Composto: considera 
o lote L como o conjunto de N grupos G, 
de fragmentos, cada um dos quais suposto 
indivisível e participando na amostragem 
com probabilidade P , bem definida, de ser 
colhido. 


41.2 — Probabilidade de colheita de um fragmento 
Os dois acontecimentos: 
* presença de F, na amostra modelo D,. obti- 
da pela reunião dos incrementos modelo; 
* toma de F,, quando F está em D,. 


são independentes, Donde, pelo teorema das proba- 
hilidades compostas: P — né * P” em que: 
1 | 


P — probabilidade de F. pertencer à amos- 


tra real: 


P' — probabilidade de F. pertencer à amos- 
tra modelo (probabilidade de «décou- 
pe»); 


P. — probabilidade de F, pertencer à amos- 
tra real, se pertencer já à amostra 
modelo (probabilidade de «tomas ). 


(1) — probabilidade de «découpe» Pi. 


Os dois acontecimentos: 
* presença de F em dx, 
* selecção de x para sede dum incremento; 


são independentes. 


Donde: 


P'= [ f(x). gíx) dx 


cj D, 
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f( x) — densidade de probabilidade de pre- 
sença de F. em x— caracteriza a 
distribuição espacial de F, e é habi- 
tualmente um dado do problema; 

g(x) — densidade de probabilidade de que x 
seja escolhido para sede dum incre- 
mento pontual — caracteriza a distri- 
buição espacial dos incrementos(1). 


A «découpe> cai no âmbito do modelo pro- 
babilista composto. 


(ii) — 'Toma. 


Pode-se admitir que F, é submetido 
individual e independentemente à operação 
de toma. Esta será pois estudada no âm- 
bito do modelo probabilista simples(?). 


4.1.3 — Modelo probabilista simples(?) 


As características Ny Ms A ea, da amostra 
E, são variáveis aleatórias. É objectivo dos capítulos 
seguintes estudar as suas distribuições, i.é., definir 
as suas leis de probabilidade e estimar os seus pri- 
meiros momentos. Para isso, suporemos que se re- 
pete independentemente a mesma amostragem um 
certo n.º K de vezes, após reincorporação. 


Seja: 


L,. — população dos fragmentos formada 
pela reunião das K amostras; 


K, — efectivo de repetição do fragmento F, 
em Ls 


P —a probabilidade de colheita de F.. 


Í 


41.4 — Distribuição de K, em Ly 
K, segue a lei binomial de momentos 


m(K)=K P, e c(K)=K P(1-P) 


Quando K — a: 


a) P. não é vizinho nem de 0, nem de 1 —a lei 
binomial tende para uma lei normal; 

b) P «1—a distribuição de K tende para log- 
-normal(*); 

ce) 1-P «1 —a distribuição de K =K-K, ten- 
de para uma distribuição log-normal. 


(1) «Découper correcta: P' =constante, 
(2) Toma correcta; P'' =1. 

(3) Vidé Bib. 11 (pgs. 337 a 365). 

(4) Vide Matheron — Bib. 15. 

(5) Desvio padrão relativo. 
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4.1.5 — Distribuição do efectivo Ny de Ep 


m(Nç) = 2 Po; cANÇ) = à P(I-P,) 
i 1 


41.6 — Distribuição da massa My de Ex 


m(M,) = ba M, P,; (Mo) ee bo M/P(1-P,) 
i i 


4.1.7 — Distribuição de Ay em Ex 


mAV=S 4 P=5 4 MP, 
i i 


S (A) = 2 A/P(I-P) - 2 aCM; (1-P;) P, 


| | 


4.1.8 — Distribulção do teor ay de Ey 


O teor a. é igual a 


Conhecem-se as leis e os momentos de A, e de 
M,- O seu quociente, no caso geral, segue uma lei 
complexa. Tentemos particularizar: 


(i) —P, próximo de 0,5: 


A e M, seguem leis praticamente nor- 
mais. O seu quociente segue uma lei que 
tende para uma lei normal quando o coe- 
ficiente de variação(*) U(M,.) tende para 
zero (M,. — cte). 


(1) — P, é vizinho de O ou de 1: 


A e M, seguem leis normais, i.é., log A, 
e log M, seguem leis normais. 

Como log a. = log A. — log M,, pode-se 
admitir que a, segue também uma lei 
log-normal. 


K 


1To 


Com a ajuda de variáveis auxiliares aditivas 
similares às que se utilizaram na análise de hetero- 
geneidade, chega-se às seguintes conclusões(!): 


a) — Em primeira aproximação, os momentos do 
teor a, da amostra E, são: 


24 “MP, 
i 
mia )=—.—=. TD = =4HA 
Ê WM P m(M,) 
“ i i k 


m( A,) 


Rr 
à (a-a,)2 M? P(I-P) 


o(a) = 


Na prática, esta primeira aproximação 
é sempre suficiente, 


b) — Em segunda aproximação haverá que adi- 
cionar algebricamente os seguintes termos 
correctivos: 


di (a-a) M2 P (1-P) 
1 
m(a,.) — — ; - 
(Om, 
i 


cr , 
25 (a-a)2 MÍP, (1-P) (1-2 P) 
ia) —> — Ea 
(2% 7 
i 


Verifica-se que os termos correctivos são sempre 
de ordem mais elevada em M /M (2) que é sempre 
um infinitamente pequeno, pelo que são desprezáveis. 


4.1.9 — Casos particulares 


4.1.9.1 — Amostragem correcta 


Por definição: y FeL:P =P, 
Donde, a partir de 4,1.8 teremos, em 1.º aproxi- 
mação: 


m(a,) =a. 


y | 
1-P à (a-a,)2 M2 
i 
o(a) =— 
P M,? 


Uí(a,) = ci(a,)/a/? 


(1) Vide Bib. 11. 


Nota-se que, quando P — P, em 1.º aproximação: 
(i) — A amostragem não é enviezada,; 


(ii) — A variância relativa UZa,) é igual à va- 
riância de erro fundamental(*): 


l-z du (a—a,)? Mt 
i 


oc EP) = — 


E a? M/? 
com a única diferença de que na 1.º, inter- 


vém as probabilidades, e na segunda as 
taxas de amostragem q. 


Vê-se pois que o modelo probabilista desenvol- 
vido neste capítulo é apenas uma generalização do 
modelo equiprovável já desenvolvido por P. Gy em 
trabalhos anteriores. 


4.1.9.2 — Fragmentos com massas iguais 


(i) — amostragem incorrecta (% F ec L; PP) 


m(aç), = bá &; /5 P&fa: 
| i 


(il) — amostragem correcta (% F. e Li: P. = P) 


m(a,.) = a, — rigorosamente justa. 


4.1.9.3 — Fragmentos com teores iguais 


m(aç) = 
ca) =0 


4.1.9.4 — Massas e teores diferentes mas com correlação 
nula 


(1) — Se a amostragem é incorrecta, é envie- 
zada; 
(il) — Se é correcta é rigorosamente justa. 


4.1.10 — Exploração prática dos resultados 


a) — Nos parágrafos precedentes considerámos 
conhecidas as distribuições de M, a e de P 
de cada um dos N, fragmentos F, que cons- 
tituem o lote L a amostrar. 

Na prática é geralmente impossível nu- 
merar todos os fragmentos dum lote e mui- 
to menos conhecer a sua massa e teor 


2 M=-m(M)=EM, P/; M,/M, — quociente entre a massa do fragmento genérico F, e a esperança da massa da 


amostra, 
(3) Vide Bib. 7 — Cap. 4.1. 
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individuais, ou ainda atrbuir a cada um 
deles uma probabilidde P de colheita 
(exceptua-se os casos seguintes: F, e L; 
P=P; P=d; P,=0). 

Assim, na prática rodeia-se esta difi- 
culdade à custa da seguinte aproximação: 
divide-se o lote em classes granulodensitá- 
rias L |. e admite-se que cada fragmento 
F. possui as características do fragmento 
médio E as da classe homogénea Lo g * que 
pertence, As características do fragmento 
médio são fáceis de determinar por pesa- 
gem e análise química duma amostra alea- 
tória de p fragmentos. 

Também se admite que todos os frag- 
mentos de L | são amostrados com a mes- 


ma probabilidade de toma Ea (1). 


bj — Considerando que este assunto já foi tra- 
tado com todo o detalhe em publicações 
anteriores de P. Gy, limitar-nos-emos a 
apresentar os resultados e unicamente no 
caso da toma ser correcta e em 1.º aproxi- 
mação: 


(1) — P'« 1 (Processo de implantação de 
incrementos): 


m(a,) = 2 dj My /My= 4, 
[ed 


o(a) = [2 ar Pias? la a pa, | fem, 


(il) — P'* = 0.5 (Processo de Partilha) 
m(a,)=4 
ta) =| E (Mud MM, | [Me 
o ia 


c) — A estimação experimental da massa M, | 
do fragmento médio E. da classe L pode 
fazer-se por: 


(i) — pesagem de M,., e contagem de N o 


(il) — pesagem e numeração de uma amos- 
tra representativa de La; 


(iii) — estimação indirecta: 


ES Ene 3 
le , Ha q, 


(nm P 
Re] 
(2) d ed 
qo a 
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é igual à esperança da taxa numérica real E 


em que: 
f — factor de forma de L (É = [ — 0.5) 
[rd EL E : 


|, — massa volúmica média dos fragmentos de 
L (ou por medida ou nor estimação em 
face do conhecimento da composição mine- 
ralógica de L) 


d. — cúmulo superior de E iu 
a =(a +a! )/2 (2). 
[4 Nu E, 


4.1.11 — Modelo Probabilista composto 


Considera-se o lote L como o conjunto dos N 
grupos G pelos quais estão repartidos os N, frag- 
mentos. 

Estes N grupos G podem ser, por exemplo, todos 
os incrementos — modelo possíveis, 

Pode admitir-se, para as leis usuais de integra- 
cão e quando o variograma da distribuição dos N 
grupos é «plat», que cada grupo G é submetido indi- 
vidual e independentemente à operação de «découpe», 
com uma probabilidade p' : 


CG ecumP =P 


Considerando que a probabilidade de colheita dos 
grupos é uniforme e que a «découpe» é correcta, 
obtêm-se os valores dos momentos e das caracterís- 
ticas da amostra E indicados nos quadros 4 e 5. 


Note-se que: 


a) — A variância UM.) não é nula porque o 
número de incrementos Q,. é uma variável 
aleatória. 


b) — Introduzindo Q = m(Q,.) — número médio 
de incrementos, obtém-se, para M =M, /N 


que é a expressão da estatística clássica 
referida no modelo de integração, o que 
significa que os dois modelos são similares 
na sua charneira: o variograma «plat», 


4.2 — Critérios de apreciação de uma amostragem na 
óptica probabilista 


Este parágrafo é uma réplica do parágrafo 3.3. 
As condições práticas são: 


(1) — Para que uma operação probabilista seja 


da classe granulo-densitária L 


são os limites superior e inferior da classe granulométrica L 
Le 
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Nd 
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dy 
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8 N vid eN vd 
ciais Ci E EEE 0 OAIJDSLIOD OULIS, 
= u E u ly 1091 
s(e-v)z (ag-r Ca-t)z pUe-s)z (az-u) (a-tz 
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E ia a À 
MN 
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78 E d-1 
«AN | ) fi 
da | N/,W (a-T)d Wa 
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exacta é necessário e suficiente que a ma- 
téria a amostrar seja dotada de homoge- 
neidade de constituição o que só é possi- 
vel idealmente; 

(ii) — Para que seja precisa (sem ser exacta), 
é necessário: 


* que seja rigorosamente justa o que im- 
plica(!) 


— que seja correcta 


YE eL;P=P 
YFEL;P=0 


— que seja satisfeita uma das seguintes 
condições suficientes: 


* MPF eL;M=M/N, 
*p (a, M)=0 
* que seja fiel; 


(ii) — Para que seja representativa (sem ser 
exacta nem precisa) é necessário: 


* que seja aproximadamente justa, L.êé., 
correcta; 
* que seja fiel. 


4.3 — Erro fundamental e erro de segregação 


4.3.1 — Introdução 


No âmbito do modelo de integração a desordem 
local apenas intervem na sua forma global através 
o termo V.,, do variograma Via (0). 

Este termo, de conteúdo probabilista, vai servir 
de charneira entre o modelo de integração regulari- 
zado e o modelo probabilista composto, em cuja ópti- 
ca nos colocâmos. 

Como, na prática, os incrementos nunca são mo- 
nofragmentários, é lógico considerar os grupos G de 
fragmentos disjuntos e complementares em L, no seio 
dos quais cada fragmento E; está indissoluvelmente 
ligado a uma estrutura fixa — estamos no âmbito do 
modelo probabilista composto. 

A integração consistirá em seleccionar, segundo 
uma lei, um certo número Q destes grupos indisso- 
ciáveis (P=Q/N). 

Ver-se-à ao longo deste capítulo, como aparece 
a integração na perspectiva descontinua da matéria. 


4.3.2 — Momentos do erro de integração El, na perspec- 
tiva probabilista 


a) — Admita-se que a(x)=a, +a.,(x), função do 
tipo da que se observa na retoma dum par- 
que de homogeneização. 


A ausência de regionalização implica 
que as populações dos teores a(x ) e a, 
são normais ou que não existe correlação 
entre o teor a e n. Neste caso as três leis 
de integração são equivalentes entre si e a 
uma tiragem ao acaso duma população 
qualquer. 

Podem pois aplicar-se os resultados do 
Modelo Probabilista à toma de Q grupos 
da população-mãe dos N grupos constitu- 
tivos de L, 

O teor a, é uma variável aleatória 
cujos momentos foram estudados no Capi- 
tulo 4.1. Os momentos de El, não são mais 
do que os momentos centrados relativos 
de Ap 

Na prática, para as três leis usuais de 
integração e sempre que D,/0=Q inteiro 
teremos; 


1. modelo de integração (vide Cap. 3). 


m(EI,) =0 
o (El,) = Vu,/Qa,2=c2(a,)/Q ay? 


2. Modelo probabilista. 


m (EI) = O 


“ 


| 1-P | 
(E) =—— dy (Ba)? M,2/aç? My? 


A 1." aproximação é suficiente (2). 


Note-se que quando P« 1, M=M /N e 
P-Q/N, teremos c2(EL)=c?(a )/Q a? 
idêntica à do modelo de integração. 


bj) — Explicitação da expressão da variância de 
EL, 
A variância de El, é igual à variância 


de “E, heterogeneidade da amostra E. 


11. N 


2 +n) (3 
Q N NA” 


Como o AE) = ( 


(à 1 No. pe 
“AO Ens 24 MA) 
Hã 


e por outro lado H (G ) = > N AN, 
n 


Cy A Justica rigorosa pode ser obtida quando as probabilidades P, não são uniformes mas são independentes da 


personalidade dos fragmentos: p(A,, P/)=0 e g(M,, P)=0. 
(2) À 2.» aproximação está indicada no quadro 4. 
(3) Da estatística clássica. 
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Se y G CL>N=N/N 


teremos H,(G )= 2 A/N = e2(A,) 


o que tomando em consideração o que se 
disse em heterogeneidade local conduz a 


nr a—a, 2) 


em que y e E têm os significados habituais. 


Hot ad) (Md 


Pode pois concluir-se: 


Considerando que Q/N-P e que N »1 


G+ro 3 | E = |. 


o? (E1,) = 
do é 


em que: 


asa tá =N/Q-1 
Se a 


Factor de amostragem que apenas depende 
de P (suposto uniforme) ou de +. A dimen- 
são e o número absoluto dos incrementos 
não intervêm. 


* 1+yE Factor de segregação e de agrupa- 
mento; 
E=0 — segregação nula; 
yY=0 — quando o efectivo de G € 1. 


“* s a, — dr, M, Fo — Factor de com- 
a M 


É posição da maté- 
ria a amostrar. Está ligado à heteroge- 
neidade de constituição e é igual a 
H. (L)/N,2. 

Como se vê, c2(EI,), na prática, nunca se 

anula, mas pode diminuir-se quando: 


(1) — Se colhe uma amostra o maior possi- 
vel; 

(ii) — se homogeniza o lote antes da amos- 
tragem; 

(dii) — se efectua um grande número de 
pequenos incrementos em vez de um 
menor número de grandes incremen- 
tos; 

(iv) — fragmentando, na medida do possí- 
vel, o lote a amostrar. 


4.3.3 — Erro Fundamental 


4.3.3.1 — Definição 


Designa-se por Erro Fundamental o valor míni- 


mo de El, quando a amostragem é correcta e y £=0. 
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El, reduz-se a EF em dois casos. 


(1) — £=0 — Distribuição naturalmente homogé- 
nea. Os incrementos podem ter 
qualquer dimensão-amostragem 
quase-eg uiprovável : 

(ii) — y=0 — incrementos de fragmentos isola- 
dos. A sua distribuição pode ser 
qualquer-amostragem equiprovável. 


Estas condições correspondem às do modelo pro- 
babilista já estudado. 


4.3.3.2 —— Momentos de EF 
1.º aproximação: 


M(EF),=0 


c(EF),= 


— Eu, | 


Z* aproximação: 


M(EF),=M(EF), 


Ae]! 
dr M 

E 

c(EF),=0(EF), 


2(1-P) (1-2 P) fe 2 /M4 
Rá P2 2 OM M 


Note-se: 


(1) — O erro fundamental é um erro patamar 
função da heterogeneidade de constituição 
da matéria e que constitui um mínimo 
incompressível. De todos os erros de amos- 
tragem é o único que nunca se anula, É o 
erro que se comete quando a amostragem 
é realizada em condições ideais. 

(ii) — O viés fundamental M(EF), é desprezável 
em face do desvio padrão 5?(EF') e anula- 
-se quando: 


— P=1 
— 8,=a, 


—M,=M,/N, 


4.3.3.3. — Condições de minimização de cEF) 


São as condições (i) e (iv) de 43.2, i.é., aumen- 
tando P ou fragmentando o lote. 
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4.3.4 — Erro de Segregação 


4.3.4.1 — Definição 


É a grandeza ES cujos dois primeiros momentos 
são iguais à diferença entre os momentos correspon- 
dentes de EI, e de EF. 


4.3.4.2 — Momentos de ES 


Indica-se apenas a 1." aproximação que, na prá- 
tica, é sempre suficiente: 


m(ES),=m(EL),-m(EF) =0 

SH ES),=cEL),-cHEF),=y £ q MEF), 

c(ES) só pode ser estimado a partir de 
c2(EL,) — análise variográfica — e de c2(EF), cuja 
estimação «à priori» e «posteriori» iremos ver na 
secção 4.3.5. 


4.3.5 — Estimação prática dos momentos do erro funda- 
mental 


4.3.5.1 — Existem fundamentalmente três métodos para, 
na prática, estimar a variância do erro funda- 
mental: 


(ij) — Um método preciso baseado na análise 
granulodensitária completa detalhada do 
material a amostrar. Por já estar sufi- 
cientemente descrito na bibliografia e por- 
que exige uma experimentação longa, de- 
licada e muito cara não o descreveremos 
em detalhe. 

(ii) — Um método rápido derivado do precedente. 
Permite fixar em poucos minutos a ordem 
de grandeza da variância fundamental 
sem qualquer suporte experimental. Apli- 
ca-se unicamente aos minérios contendo 
um elemento crítico sob a forma dum mi- 
neral único. 

Pode adaptar-se a outros materiais 
a granel desde que os seus constituintes: 


* Sejam nitidamente expressos; 

* tenham composição e propriedades 
aproximadamente constantes; 

* estejam justapostos no seio dos 
fragmentos e não em solução uns 
nos outros. 

(iii) — Um método prático baseado na colheita 
de uma amostra representativa de cerca 
de 50 fragmentos da classe granulomeétri- 
ca mais grosseira e na determinação da 
massa, do volume e do teor de cada frag- 
mento. 

É um método mais simples, mais rá- 

pido e mais barato que o primeiro. Exige 
unicamente equipamento laboratorial cor- 


rente e é aplicável a numerosos materiais 
a granel: minérios, produtos metalúrgicos, 
materiais não minerais, etc. 


4.3.5.2 — Método expedito (ii) 


E Bite a—-a Mo |; 
Seja S'=M, X o | | 
eme 


Da exploração dos resultados da aplicação da 
análise granulodenaitária, tem-se: 


s= E y Eta) Mg/m MO 
Er 


c(EF),=S' (1-P) /PM =S'(1/M -1/M/) (2) 


* Transformação de S' num produto de dois 
factores. De (1) Ltira-se 


! — E o Ea F. E 
= Do 2 VaPo (Poor)! Mg, 
x 8 


A experiência mostra que o teor a, depende 
muito mais de q do que f, pelo que se fará a 1.º 
aproximação: 

de mta e ido 


Verifica-se também que M af /M é habitualmente 
varia pouco duma classe densitária para outra o que 
conduz à 2.* aproximação: 


VLoCL:Mo/M=M o /M 


Então: 


Ss = S', X S', em que 
o — Er Fi E 
3, = 2 Ha (aç-a,) M./a, se 
Ss E 2” ME /M, 
[e 


* Expressão simplificada de S'. 
Os valores extremos de Ss", são: 


O — Quando o lote é rigorosamente homogêneo. 
c= [(1-a))/2,]*a, Ha + &, UH; quando o lote é 
constituído por particulas inteiramente livres. 
Na prática: O< S',<c. 
Escrever-se-á então: 
é = 61 
c — Parâmetro mineralógico [M L“1(1) 


1 — Parâmetro de libertação (sem dimen- 
sões). 


(1) É necessário exprimir os teores em valor decimal e em mineral e não em metal, 
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* Expressão simplificada de S', 
E 3 
S,= > E, d,M,/M, 
x 


Como di diminui rapidamente pode considerar-se 
unicamente o d da fracção mais grosseira q 
= dg = d. 


“= f.g. di em que 


f — parâmetro de forma (sem dimensões) = 0.5 
g — parâmetro granulométrico (g = 0.25) 
d— dos 


1 1 
* gHEF) =cl f.g. d É En ) 


Ma My 


c.1.f.g. d3'/M, 
Cd/M 


um 


Já suficientemente descrito na bibliografia. 


4.3.5.3 — Método prático 
Procede-se como indicado em 4.3.6.1 (iii). 


(i) — Estimação de S', 


; 2 airê Z 
Ss, — > (a—&o) M,/Va do Mo 
q 


em que 


a —teor de F 
q q 


M, — massa de E 
Na — volume Fr 
ão — teor de Lo 
M, — massa de Lo 


(ii) — Estimação de S', 
E = rs E ie 
= > VM /M,=t Sai M 
a 
= 3 
ca o e M/M, 
E) 


ou então S,=f.g. d'= V/4 


o que conduz a: 


MEF) = 8, XS)/M =8! V/4M, 


(1) De notar que C é função de d a partir de 1. 
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4.3.6 — Exemplos de aplicação 
4.3.6.1 — Método expedito 


O método descrito em 4.3.5.1 permite resolver 
problemas do tipo: 


(1) — Cálculo da variância fundamental 
(EF) = € d3/M 


(ii) — Cálculo da massa da amostra a escolher 
(C e d fixados) para que qcHEF) < q,? 
M, > Mpo = €C d3/0,2 


(iii) — Cálculo da dimensão crítica a que é ne- 
cessário moer o minério nas condições do 
problema posto 


d<d =VM 0,2/C (1) 


4.3.6.2 — Método prático 


Apresenta-se um exemplo dum minério de ferro 
da região de Cassinga, Sul de Angola. 


De acordo com o método operatório já sumaria- 
mente descrito constituiu-se um lote de trabalho L, 
do qual se colheu, ao acaso, uma amostra de 46 
fragmentos da classe granulométrica mais grosseira. 

Para cada um destes fragmentos foi determi- 
nado: 


— massa M, [£] 
— volume Va [em3] 
— teor a. [SG Fl 


Os resultados constam do quadro 6. 


Seguidamente calcula-se: 


aq = 389% F,: as = 1487.17 (% F))? 


M,=5971g ; V5=1733em'; V= 37,67 cem 
2 a 
ba (a —aç)? M/Va aq Mo = 488.5 g? ema 


ds me ao 2 
S, = b3 (a-ag)” M/V, & 


o Mo = 0.63g. cm: 


S, = 0.25 V = 9,4 em: 


. 


4 
&' 


é E F é. 
S, XS,=59g » cHMEF) = 5.9 | E ado ) 
A partir destes resultados é possível calcular a 


constante de amostragem C para se poderem resolver 
os três problemas clássicos. 
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QUADRO 6 


q a q 
grama cm! Fe % grama cm 
163 42 55,31 189 49 
121 | 44 21,33 157 48 
116 39 37,86 116 
127 39 45,62 66 24 
78 29 13,91 141 53 
102 26 50,86 128 39 
60 23 42,65 138 50 
142 48 44,13 Ss 31 
105 34 26,91 150 42 
119 45 36,84 100 41 
61 19 43,22 208 60 
149 57 4,11 127 39 
95 37 6,04 146 49 
95 89 7,38 117 30 59,07 
198 48,91 102 38 20,53 
107 53,14 mn 29 10,38 
172 43 61,58 81 29 44,82 
120 29 62,72 82 36 3,40 
119 33 38,20 105 31 41,28 
104 28 44,48 116 33 48,13 
121 33 44,93 110 31 46,83 
EE] 32 47,56 67 26 7,64 
148 40 63,52 
24 52 20 7,80 5371 1733 37,91 
q 
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Assim, considerando f=0.5, teremos: 


d3 = Vy/f = 75.34 em 
d = 4.22 cm 
C=S/d3=0.078g.em: 


1 1 
c2(EF) = a a wo d3 


O cálculo do desvio padrão conduziu aos seguin- 
tes valores: 


Massa da amostra c(ap) [% Fe] E Zon 
10 kg 1,09 2,18 
100 kg 0,35 O,T0 
1000 kg 0,11 0,22 
10 000 kg 0,035 0,07 
4.3.7 — Notas 
( dg 1 ) 
Na expressão c(EF) == ww J|ef.lpg.d3 
Mp M, iá 


c — Parâmetro mineralógico — traduz a influên- 
cia da constituição da matéria amostrada. 
É devido a este parâmetro que quando se 
utiliza esta expressão se têm que exprimir 
os teores em valores decimais de mineral 
e não de metal; 

| — parâmetro de libertação traduz simultanea- 
mente a influência da estrutura da matéria 
e do grau de individualização dos seus cons- 
tituintes; 

f — parâmetro de forma — traduz a influência — 
fraca — da forma dos fragmentos; 


g — parâmetro granulométrico — traduz a in- 
fluência — fraca — da maior ou menor gra- 
nularidade do minério; 

d — da, — traduz a influência do estado absoluto 
de fragmentação. 


44 — Análise do Processo de materialização dos incre- 
mentos pontuais 


441 — Introdução 


Como vimos em 1.3.1, o processo de amostragem 
por implantação de incrementos («colheita») pode ser 
considerado como a sequência de quatro operações: 
integração, «découpes, toma e reunião. 

A integração foi estudada no Cap. 3. A reunião 
é uma operação exacta. A materialização dos incre- 
mentos pontuais, efectuada em duas etapas «décou- 
pe» e toma, é geradora de erros que são estudados 
no âmbito do modelo probabilista. Antes de entrar 
no detalhe, iremos examinar o conjunto de transfor- 
mações que permitem passar do modelo à realidade 
e quais os erros que originam, 
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44.2 — Lotes com uma dimensão temporal (OT,) 


(ij) — Ponto de partida — incremento pontual. 
O incremento é representado por um 
ponto x do eixo dos tempos. O teor a, 

de L é rigorosamente: 


FT ft 
& = [ ate) lx) ox / | oleo dx 
"0 


Lo | 


(lj — Primeira transformação — incremento 
segmentado. 

O incremento é representado por um 
segmento de comprimento T centrado 
em x. 

O teor do material amostrado é repre- 
sentado pela função a (x) regularizada 
de a(x) à escala do segmento de compri- 
mento T. E a média de a(x), ponderada 
pelo caudal pontual u(x) no interor de T.. 

O caudal é representado por ju (x), 
regularizado de u(x) à escala de T. E a 
média simples no interior de T.. 


O teor a, do lote será: 


Wi T, | Tt, 

m 

a =, faico Wu. (x) dx [ultimo dx 
E Õ E] LO o 


como T.«'T, pode admitir-se 
T, T, 
a = fait (x) dx fico dx 
“0 o 


Esta identidade prática não é gera- 
dora de erros apreciáveis, desde que 
Ev agi 


(iii) — 2." transformação. 
A matéria, que suposemos condensada no 
eixo dos tempos, ocupa realmente um cer- 
to volume no espaço geométrico contido 
entre dois planos paralelos (mas não pa- 
ralelo: ao eixo dos tempos) distanciados 
de T, 

É a «découpe» modelo, fatia de espes- 
sura constante da matéria não discreti- 
zada. 

Esta identidade não é geradora de 
eTTOS. 

(iv) — 3.º transformação — «Découpe» real — 
Erro de «Découpes». 

A operação de colheita é habitualmen- 
te realizada por um recipiente —o amos- 
trador — que no instante da colheita é 
levado a interceptar a envolvente do fluxo 
de material recortando naquele fluxo um 
volume a que chamaremos «découpe» real. 
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A experiência mostra que nem sempre 
as superfícies são planas nem paralelas 
o que conduz a uma não uniformidade de 
probabilidade de colheita dos diferentes 
fragmentos, devido à maior ou menor 
correlação entre as características granu- 
lo-densitárias de um fragmento e a sua 
posição espacial no fluxo do material. 

É esta incorrecção que gera os erros 
de «découpe». 


(v) — 4," transformação — «découpe» discreta. 


A «découpe» discreta é constituída 
pelo conjunto de fragmentos cujo centro 
de gravidade está no interior da «décou- 
pe» real. Corresponde à selecção que se 
opera, em condições ideais, quando os 
fragmentos caem sobre uma aresta fixa 
ressaltando para o lado onde se projecta 
o seu centro de gravidade. 

Esta transformação corresponde a 
uma equivalência estatística, pois pode 
demonstrar-se que a diferença entre as 
composições das «découpes» real e dis- 
creta é uma variável aleatória de média 
nula. 

Não é portanto geradora de erros. 


(vi) — 5.º transformação — incremento real — 
erro de toma. 

O incremento real é constituído pelos 
fragmentos que penetram efectivamente 
no amostrador e nele permanecem. 

Às condições que conduzem a uma 

s identidade entre a «découpe» discreta e o 
incremento real são: 


a) — Condições ligadas à queda dos frag- 
mentos: 
* os fragmentos devem cair um a 
um; 


* devem cair verticalmente; 


* devem cair sem energia de «spin+. 


b) — Condições ligadas ao amostrador: 

* a abertura do amostrador deve ser 

limitada por arestas lineares sem 
espessura; 


* as arestas devem ser horizontais; 
* as arestas devem estar suficiente- 


mente afastadas uma da outra, 
sem o estarem infinitamente; 


* as arestas devem estar imóveis; 


* o amostrador deve ser concebido 


por forma a respeitar a integri- 
dade dos incrementos. 


Quando uma destas condições não for 
respeitada, não se verifica a regra do 
ressalto, não havendo identidade entre a 
e«découpe» discreta e o incremento real, 
gerando-se assim um erro de toma. 

Pode pois concluir-se que na materia- 
lização dos incrementos reais há duas 
operações susceptíveis de perar erros: 
«Découpe» e Toma. 

A Pig. 4 esquematiza o que acabamos 
de escrever. 


RECAPITULAÇÃO ESQUEMÁTICA DA PASSAGEM 
DO MODELO DE INTEGRAÇÃO A REALIDADE 


Projecção a duas dimensões 
(Segundo P. Gy) 


“a 


M 


| - incremento pontual 


———————T——. 


2 = Incremento segmentado 


- Decoupe" modelo 


4 - “Decoupe" real 


5 - “Decoupe" discreta 


000065 
000. 


& - Incremento real 


Fig. 4 
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MARK 1 TIME-SHARING SYSTEM 


Para além da Prestação de Serviços às Empresas nossas clientes, 

A colaboração e apoio que, ao longo de oito anos, temos prestado às Universidades, 
Institutos Técnicos, Liceus e Centros de Investigação, constitui uma linha de força 
maior da nossa actividade. 

Tal colaboração tem-se traduzido, não só em palestras, cursos de análise e pro- 
gramação, estágios, resolução de problemas académicos, etc., como em condições 


especiais de utilização dos sistemas Time—Sharing. 


GE — 265 MARK I 
(Sociedade Time Sharing) 


DEC SYSTEM 10 
(Lab. Nacional de Engenharia Civil) 


A instalação de terminais nos locais de trabalho, permite o contacto telefónico com 


os dois computadores acima referidos e, em casos especiais, com a rede europeia de 


MARKTII 
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44.3 — Lotes com 1, 2 ou 3 dimensões espaciais 


Pode reconstruir-se, para cada um destes lotes, 
a sequência que acabámos de indicar. 

Pode ainda acrescentar-se que, suprimindo a 4.º 
transformação («découpe» discreta) e modificando 
ligeiramente a formulação da 5.º, podem deduzir-se 
um certo número de conclusões válidas para sólidos 
compactos não fragmentados tais como 


* jazigos mineiros; 
* produtos metalúrgicos: lingotes ou chapas 
metálicas. 


4.5 — Erro de «découpe- 
4.5.1 — Definição 


Como vimos, o erro de «découpe» surge quando 
a «découpe» real, efectuada pelo amostrador no do- 
mínio D, ocupado pelo lote L, difere da «découpe» 
modelo definida em função do número p de dimen- 
sões do modelo representativo do lote. 

Em coordenadas cartesianas: 


* p=3— é um cubo; 

* p=2— é uma coluna de secção constante 
e quadrada e de altura igual à es- 
pessura do lote; 

* p=1l — é uma fatia de espessura constante 
cortando inteiramente a seccão do 
lote. 


4.5.2 — Condições teóricas de amostragem correcta(") 

Uma amostragem é correcta quando assegura a 
todos os fragmentos F. de L uma igual probabilidade 
de pertencerem à amostra E. 


Ora, a probabilidade de colheita de É; é: 


P, = fixos g(x) dx 
D, 


em que 


f(x) — densidade de probabilidade de presença 
de F, no ponto x — distribuição espacial 
de F.; 


g(x) — densidade de probabilidade de localizar 


um incremento no ponto x(2) — distri- 
buição espacial dos incrementos. 


A condição de correcção será pois: 


WWF eL: P= faco gix) dx=P 
o“D 


(0<P<1) 


L 


As condições teóricas de distribuição dos incre- 
mentos g(x) que asseguram a correcção, para os 
diferentes tipos de distribuição espacial, ou seja, de 
homogeneização do lote, estão expressos no quadro 7. 


QUADRO 7 


DISTRIBUIÇÃO CORRECTA DOS INCREMENTOS 


Distribuição Distribuição dos incrementos, 
Tipo de homogeneidade dos fragmentos assegurando amostragem Tipo de Lote 
F,(x) correcta 

Tridimensional | f, g(x, y, Z) — qualquer de forma qualquer 

bidimensional f(z) g(x, y) — uniforme em z tabular 

unidimensional f(y,Z) g(x) — uniforme em (y,z) | parqueamento alongado; 
lote em deslocamento 
(correia transportadora) 

de revolução E (pz) g(g) — uniforme em p e z | de revolução 
(enquartação ) 

heterogéneo f(x,y,2Z) £., — uniforme de forma qualquer 


() Por amostragem considerámos aqui o conjunto «integração» — «découpes, 
(2) A probabilidade de que D, englobe um ponto qualquer de D, é DD, e é igual a 5 f D g(ix) dx. 
E L. L 
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Pode concluir-se dizendo que: 


* Se podem seleccionar livremente a(s) 
coordenada(s) de homogeneidade; 


* se deve seleccionar a(s) coordenada(s) de 
heterogeneidade e atribuir-lhe(s) uma den- 
sidade uniforme de probabilidade. 


4.5.3 — Condições Práticas de «découpe» correcta 
4.5.3.1 — Lotes com uma dimensão temporal 


4.5.3.1.1 —= Amostradores «transversais» 


Definem-se amostradores «transversais» como 
aqueles que colhem a totalidade do fluxo durante uma 
parte do tempo. 


Demonstra-se que um amostrador «transversal» 
realiza uma «découpe» correcta quando e só quando 
assegurar a cada linha de fluxo y de material, uma 
igual duração de colheita Tl Y) = constante. Caso 
contrário a «découpe» é incorrecta. 


Também se demonstra que o erro de «découpe» é 
identicamente nulo quando e só quando a «découpe» 
realizada é correcta, 

As condições práticas da correcção estão es- 
quematizadas no quadro 8. 

O respeito pelas condições geométricas do amos- 
trador está esquematizado nas Fig. 6 e 7. 


4.5.3.1.2 —= Amostradores «não transversais» 


São amostradores que colhem uma parte do fluxo 
durante a totalidade ou uma parte do tempo. 

Originam erros sistemáticos que podem ser gra- 
ves. São desaconselháveis e estão a cair em desuso. 


4.5.3.2 — Lotes com dimensão zero 


Um lote tem dimensão zero quando está natural- 
mente dividido num grande número de unidades de 
peso estatístico igual e de peso físico aproximado 
(ex.: comboio de vagons de minério — método da 
pazada, etc.). 

O erro de «découpe>» é nulo se se proceder à 
amostragm primária desses lotes através a colheita 
dum certo número de unidades inteiras. 


4.5.3.3 — Lotes com uma dimensão espacial 


Só pode ser considerada correcta a «découpe» 
do tipo da que é executada com um quadro sobre 
uma correia transportadora parada (vide Fig. 5). 


4.5.3.4 — Lotes de material a granel com 2 ou 3 dimen- 
sões espaciais de revolução 


É extremamente difícil na prática efectuar uma 
e«découpe» correcta por muito onerosa ou mesmo im- 


possível, 


4.5.3.5 —- Conclusão 


Pode concluir-se que só a amostragem dos lotes 
em deslocamento continuo (unidimensional) ou des- 
contínuo (dimensão zero) pode assegurar, na prática, 
a realização de uma «découpe» correcta. 


4.6 — Erro de toma 
4.6.1 — Definição 


O erro de toma surge quando o conteúdo do 
incremento real difere do da «découpe» discreta, isto 
é, quando a regra do ressalto não é respeitada, 

Vejamos seguidamente as condições práticas de 
toma correcta (1). 


Fig. 5 


(1) As condições teóricas estão descritas em 4.4.1 (vi). 
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QUADRO 8 


AMOSTRADORES TRANSVERSAIS — CONDIÇÕES DE «DECOUPE- CORRECTA 


Condição de «découpe» correcta: todas as linhas de fluxo y devem ser interseptadas num intervalo 
de tempos constante 


Uniformidade de T (y) — duração de colheita — no 


espaço a (dado 


E e— 


Uniformidade de To Y) no tempo 


= = 


Condição cinemática: 
Velocidade NAS) uni- 
forme em cada incre- 
mento 


Condição geométrica: 
Largura W (y) 

' Forma e instalação 
do amostrador 


Condição cinemática: 


Condição geométrica: 
Largura uniforme de 
um incremento para 
outro | 


Velocidade Va83) uni- 
forme de um incre- 
| mento para outro 


a) Comando eléctrico; 
Correcto se o mo- 
tor for sobredi- 
mensionado e as 
posições de repouso 
do amostrador fo- 
rem suficientemen- 
te afastadas do flu- 
xo do material, 


b) Comando hidráu- 
lico, pneumático ou 
manual; sempre in- 
correcto (a evitar). 


Concepção do amos- 

trador: 

1) Trajectória rectili- 
nea: 
correcta se arestas 
paralelas, 

2) Traject. circular: 
correcta se arestas 

radiais. 

3) Qg. outra traject.: 


(amostra. manual) | 


sempre incorrecta. 


2) Comando eléctrico: 
Praticamente cor- 
recta se as flutua- 
ções de tensão e 
frequência da rede 


eléctrica não forem | 


elevadas. Erro si- 
milar ao da flutua- 
cão do caudal ED. 


b) Comandos  hidrãu- 
lico, pneumático ou 
manual: incorrec- 
tas — a evitar. 


Largura uniforme: 


Correcta se se limpa 
periódica e frequente- 
mente a abertura do 
amostrador e se se vi- 
gia a deformação e o | 
desgaste. 


| Instalação e manuten- 
ção do amostrador: 


a) Interseptar totalmente o fluxo do material. b) Evitar as poeiras. c) Posições convenientemente afasta- 
das do fluxo do material. d) Limpar periódica e frequentemente a abertura do amostrador. e) Vigiar a 
deformação e o desgaste da abertura, 


Pode pois concluir-se que a realidade 
é estatisticamente equivalente ao modelo 
correcto. 


4.6.2 — Lotes de material a granel com uma dimensão 
temporal amostrados com um amostrador «trans- 
versal» 


(ij) — Os fragmentos tombam um a um. 

Na prática os fragmentos tombam 
sob a forma dum fluxo mais ou menos 
compacto. Quando um fragmento ressalta 
sobre uma aresta interfere com a trajec- 
tória dos outros fragmentos. 

No entanto pode admitir-se que na 
parte central do fluxo a reacção é esta- 
tisticamente simétrica, 

Na periferia a reacção é centrífuga. 
Mas como o amostrador possui duas ares- 
tas pode admitir-se que elas desempenham 
papéis simétricos. 


(ii) — Os fragmentos caiem verticalmente. 
Também aqui a alteração é simétrica 
ao nível das duas arestas do amostrador, 
podendo concluir-se que a realidade é es- 
tatisticamente equivalente ao modelo cor- 
recto. 


(iii) — Os fragmentos caem sem energia de 
«spin», 

Também aqui a alteração é simétrica 
em relação às duas arestas podendo tam- 
bém concluir-se que não há erro de toma 
nos casos reais. 
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«DÉCOUPE- CORRECTA: Respeito pela condição geométrica — Concepção do amostrador 


(Segundo P. Gy) 


1) Trajectória rectilinia 


CORRECTA INCORRECTA 


2) Trajectória circular 


o Eixo O Eixo 
nº 0 
1 
! h 
! ! 
CORRECTA INCORRECTA 


Fig. 6 


(iv) — A abertura do amostrador é limitada por 
arestas lineares. 

Na medida em que as pequenas defor- 
mações devidas ao desgaste não compro- 
meterem nem a forma nem as dimensões 
originais do amostrador, pode afirmar-se 
que não hã erro de toma. 


(v) — As arestas do amostrador são horizon- 
tais. 


Sempre que não forem horizontais a 
regra do ressalto não pode ser respeitada, 
cometendo-se um erro de toma que afecta 
tanto mais os fragmentos quanto mais 
grosseiros forem (vide fig. 9). 


(vi) — As arestas devem ser largamente afasta- 
das uma da outra, 


Se esta condição se não der, as par- 
tículas, sobretudo as mais grosseiras, po- 
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INCORRECTA INCORRECTA 


(depósito (divisor com 


sobre as um único la 


arestas) do) 

o) Eixo Eixo 

[à] à) 
INCORRECTA INCORRECTA 


dem ressaltar para fora do amostrador 
como se indica na Fig. 8. 

Para que a toma seja correcta é ne- 
cessário que a largura W da abertura do 
amostrador, no sentido do seu desloca- 
mento, seja 


vd>smm: W>W-Sd 
Yd<3mm: W>W =10mm 


Se esta condição não for satisfeita, 
ocorre um efeito estatístico de crivagem 
que altera selectivamente a probabilidade 
de toma dos fragmentos de maior calibre. 


(vii) — As arestas do amostrador são imóveis. 


Como é óbvio o amostrador tem de 
deslocar-se em relação ao fluxo do ma- 
terial. 
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«“DÉCOUPE- CORRECTA: Respeito pela condição 
geométrica — Instalação do amostrador 


(Segundo P. Gy) 


A, B, A  B 


= O  ——————— o o oO o 


INCORRECTA 


D, 


o 
e ia 


me 


a a 
Sã q io a e e 


INCORRECTA 
dr — Secção horizontal do fluxo do material 


&, — Secção horizontal do fluxo de poeiras 


Fig. 7 
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Direcção do deslocamento 


ca = 
do amostrador 


| 


5, 


Fig. 8 


A experiência prova que a toma é 
correcta se for respeitada e seguinte a 
condição cinemática: 

É necessário que a velocidade V dum 
ponto qualquer da abertura do amostra- 
dor satisfaça à relação 


V< Va W/W, 


em que 


V, =400 mm. s-1; 
W,- largura minima definida em (vi); 
W -—largura da abertura do amostrador. 


Ex.: Consideremos em minério com 
d=25mm.W=3d=7Tima. 


Para W=Wo= 10 mm — Vo= Vo 
=400 mm. s-1 


Para W=150 mm as dd Pe à fa 
= 800 mm. s-1 


A massa do incremento, sendo propor- 
cional a W/V, é a mesma em ambos os 
casos. 


(viii) — O amostrador deve ser concebido por 


forma a respeitar a integridade dos incre- 
mentos. 
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ERRO DE TOMA — INFLUÊNCIA DA ORIENTAÇÃO DAS ARESTAS 
Vistas num plano vertical perpendicular ao eixo de deslocação do amostrador 
(Segundo P. Gy) 


11 Amostrador de colher: 


A aresta AD é horizontal 


2) Amyustrador de godet; 
A aresta AD é horizontal. 


Pa 


3) Amostrador de bolso: 
4 aresta AD é obliqua. 


Há risco de alteração selectiva da probabilidade de toma dos fragmentos 
de maior calibre. 


nem a o? 
ma Ts, 


pe] 


a 


E Seb bat 


4d Amostrador de colher envolvente: 


A aresta AD é curvilinea. 
Mesmo defeito que em 3. 


5) Amostrador de pulpa de colher; 
A aresta AD é vertical. 
Mezmo defeito que em 3 mas mais acentuado. 
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Deve assegurar: 


* que um fragmento que penetre no 
amostrador não possa sair (colher com 
pelo menos o dobro da capacidade do 
incremento — profundidade > 3d; ares- 
tas com inclinação > 45º ou 60º con- 
soante o material for seco ou húmido; 
evitar correntes de ar desempoeirado- 
ras); 


* que não haja contaminação por mate- 
rial estranho (evitar fontes de poeiras, 
aço inoxidável, protecção de abertura 
com uma tampa). 


4.6.3 — Lotes de material com uma dimensão temporal 
amostrados com um amostrador não «transversal» 


É de evitar. 


4.6.4 — Lotes com 0, 1, 2 e 3 dimensões espaciais 


(i) — O erro de toma é nulo quando se procede 
à amostragem dum lote com zero dimen- 
sões através da colheita dum certo nú- 
mero de unidades inteiras, 


(ii) — O erro de toma é nulo na amostragem de 
lotes com uma dimensão espacial quando 
se procede à colheita, com o auxilio de 
um quadro, sobre uma correia transpor- 
tadora parada, desde que se tomem pre- 
cauções quanto aos fragmentos de maior 
calibre colocados sob as paredes que limi- 
tam o quadro. 

(iii) — A toma, nos lotes com 2 e 3 dimensões 
espaciais e de revolução, é sempre deli- 


cada, como já se disse ao falar na «dé- 
coupe». 


De notar que a correcção da toma é condição 
necessária de correcção da amostragem. Só lotes do- 
tados de homogeneidade de constituição (caso ideal 
sem existência real) podem ser amostrados correc- 
tamente com tomas incorrectas, 


5. PROCESSOS DE PARTILHA (DIVISÃO) 


5.1 — Definição 


O processo de partilha consiste na sequência das 
quatro operações já descritas em 1.3.2, a saber: 


* Divisão 
* Separação — conduzindo às fraccões-reais; 

* Escolha  — conduzindo aos incrementos reais; 
* Reunião -— conduzindo à amostra real. 


— conduzindo às fracções-modelo; 
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Esta sequência pode ou não ser repetitiva. 

É a escolha que constitui a amostragem propria- 
mente dita. As operações preliminares de divisão e 
separação limitam-se a transformar L num lote com 
zero dimensões, antes da sua amostragem. 

E um processo que se aplica exclusivamente aos 
lotes manipuláveis, isto é, aos lotes que se aceita 
manipular com o objectivo de amostragem ou aos 
lotes cuja transferência impõe a sua manipulação 
por razões estranhas à amostragem, 


5.2 — Momentos do erro cometidos na escolha 


Se o lote a amostrar é pequeno e a divisão e a 
separação forem correctas, pode-se considerar que o 
erro total de amostragem é igual ao erro de integra- 
ção El, cujos momentos são: 


m (El,)=0 =» m(a,)=a, 


5 (EL) — (5 = RI (a,—%) /Na, 
U(e= (qm) 06 


Para Q=N/25 U'(a,)=U'(a )/N 
em que: 


N — número de fracções reais G, cuja reunião 
constitui L,; 


Q — número de fracções seleccionadas. 


5.3 — Critérios de apreciação de uma partiha 


(1) — Para que uma partilha seja rigorosamente 
justa é necessário e suficiente que: 


* a divisão seja correcta; 

* a separação seja correcta; 

* as massas M º das fracções reais G, 
sejam todas iguais. 


(li) — Para que uma partilha seja aproximada- 
mente justa é necessário e suficiente que: 


* a divisão e separação sejam correc- 
tas; 

* as massas M sejam da mesma or- 
dem de grandeza. 


(iii) — Para que uma partilha conduza a disper- 
são nula é necessário e suficiente que 
c*(EI,) = o, o que só sucede em dois ca- 
sos limites: 


* Q=N — todo o lote é amostrado; 
* a =a, — homogeneidade de consti- 
tuição. 
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(iv) — Nem sempre é possivel realizar uma par- 
tilha fiel (1). Em certos casos seria neces- 
sário quer: 


* conservar todo o lote; 


* fragmentar e homogeneizar o lote 
antes da divisão para reduzir UZ(a ). 


* rever o nível de fidelidade R 2 


(v) — Diz-se que uma amostragem comercial é 
equitativa quando o valor comercial do 
lote, calculado a partir do teor da amos- 
tra, admite, para média, o valor calculado 
a partir do teor do lote. 

Quando o valor de L é uma funcão 
linear do teor, pode dizer-se que a amos- 
tragem é equitativa desde que o desvio 
gistemático, se existir, tenha uma igual 
probabilidade de favorecer o comprador 
e o vendedor, ou, na hipótese de uma 
amostragem de rotina repetida numerosas 
vezes ao longo do tempo, a esperança ma- 
temática do viés, considerado no tempo, 
seja nula, 


Pode resumir-se dizendo: 

Toda a partilha se pode dividir numa operação 
material de separação e numa operação imaterial de 
escolha. Estas operações só são dissociáveis quando 
a partilha dá origem a duas ou n fracções gémeas. 
A justiça é uma propriedade da separação e a equi- 
tabilidade é uma propriedade da escolha. 

Para que uma partilha seja equitativa é necessá- 
ria uma das duas condições seguintes: 


* ou que a partilha seja justa; 


* ou que, sendo a partilha enviezada, 
dê origem a duas ou n fracções gé- 
meas, entre as quais se selecciona, 
com uma probabilidade uniforme de 
escolha, a amostra, utilizando por 
exemplo, o método da tiragem à 
sorte. 


5.4 — Breve análise critica dos métodos e aparelhos de 
partilha 


5.4.1 — Pazada fraccionada (vidé fig. 10) 
A pazada fraccionada de ordem p é uma gene- 


ralização da pazada alternada, a qual consiste em re- 
tomar um lote de material com uma pá (manual ou 


(1) Partilha em 3: - fa,-4,)/Nº nº, < RR. 


Partilha em p: UlMa)/Q=E (a,—-0,)/Q Na! < Rç2 


mecânica) e colocar a primeira pazada em E, a se- 
gunda em E,, a terceira em E, e assim sucessiva- 
mente. Obtêm-se duas fracções gémeas E, e E.. 

A divisão é correcta se a taxa de enchimento da 
pá for sensivelmente constante. 

A separação é correcta se o operador proceder 
«cegamente», seguindo regularmente o contorno de 1, 
até ao esgotamento. 

A partilha é equitativa desde que coresponda 
efectivamente a uma separação em p amostras po- 
tenciais gémeas e que a escolha da amostra real seja 
feita ao acaso dentre aquelas. 

Para que a pazada seja fiel, poder-se-á calcular 
Q(quando Q « N) como segue (a partir da variân- 
cia ) 

Q > Q,=UZa,) 2 


Conhecendo aproximadamente a massa do lote 


M, e a massa média de cada pazada M., obtém-se: 
N=M,/M e p=N/Q 


Na prática e como s(EF) só depende da massa 
da amostra(p), para reduzir o erro de segregação 
(devido a y) para valores aceitáveis e compatíveis 
com a economia do processo, admite-se que cada 
fracção gêmea seja constituída por 50 pazadas. Assim, 
se p for a taxa de redução escolhida e M a massa 
do lote, a dimensão média da pazada será: 


M =M, /50p 


5.4.2 — Enquartação (fig. 11) 


A correcção da separação é difícil. A equitabili- 
dade só é assegurada quando os dois sectores circula- 
res opostos, retidos como amostra, sejam escolhidos 
ao acaso, após a separação ter sido efectuada. 


Pode concluir-se que a enquartação, ainda muito 
usada, é inutilmente complicada e cara pelo número 
considerável de operadores que exige(2). Em nume- 
rosos casos, é aconselhável a sua substituição pelo 
método da pazada. 


5.4.3 — Divisor Jones 


É o amostrador mais utilizado nos laboratórios 
para lotes de sólidos de peso até cerca de uma cen- 
tena de quilos constituídos por fragmentos com 
d<Z2 adgimm(Fig. 12). 


De notar que a, é um dado do problema e Q e N são dados de construção do amostrador. 
(2) Fidelidade: estamos no caso em que N=4 e Q=2. Teoricamente é fiel se E (a,—a,)2/16 a,2 < R,2 Vê-se pois que 


a enquartação não oferece qualquer possibilidade de resolver o problema, excepto se se aumentar o número de andares. 
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PAZADA ALTERNADA 
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il 
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1) Formação do cone -"homogenização' 


2) Formação do tronco de cone 


h 


3) Divisão em sectores 


Fig. 1 


Para que a partilha seja equitativa é necessário 
e suficiente que se proceda a uma tiragem à sorte 
para escolha da fracção — amostra real, 

Quanto à fidelidade pode dizer-se: 


(i)—o erro fundamental só depende da massa 
da amostra e, como tal, é independente 
das características do amostrador; 


(ij —o erro de segregação depende 


* da heterogeneidade de distribuição 
caracterizada pelo parâmetro de se- 
gregação £. Reduz-se £ homogenei- 
zando o lote; 

* da dimensão dos incrementos, aqui 


da fracção do lote que atravessa 
cada uma das ranhuras, 


Pode pois concluir-se que um amostrador Jones 
é tanto mais fiel quanto maior for o número de 
ranhuras que tiver. 

Dados a generalidade e o interesse da sua utili- 
zação, indicam-se seguidamente as recomendações ne- 
cessárias para garantir a correcção da divisão e da 
separação: 

1. O operador deverá despejar a pá no centro 
das ranhuras, onde estas têm a mesma pro- 
fundidade. 

2. O operador deve despejar a pá num movi- 
mento lento e contínuo para evitar o risco de 
obstrução dos canais e de emissão de poeiras. 

3. A largura das ranhuras deverá ser > 2d+5 
mm, sendo d o maior calibre das partículas 
do lote, para impedir a obstrução dos canais 
e permitir a sua limpeza cuidadosa, 


Fig. 12 
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4. Convém que o material do amostrador seja 
de aço inoxidável. Os canais devem ter forte 
inclinação (> 45º) e devem evitar-se os ân- 
gulos vivos. 


5. O número de canais deve ser obrigatoriamen- 
te par e todos devem ter a mesma largura, 
por forma a que a probabilidade para que um 
fragmento se encontre numa das duas frac- 
ções seja igual a 0.5. 


6. A experiência mostra que, qualquer que seja 
o erro sistemático cometido em cada andar, 
quando se procede a uma divisão em cascata, 
uma amostra obtida através um número 
igual de incrementos à esquerda e à direita 
do divisor tem elevada probabilidade de ser 
justa. 


5.4.4 —- Divisão em MN utilizando divisores rotativos 


(Fig. 13) 


Estes amostradores asseguram sempre uma divi- 
são e separação correctas. A equitabilidade é asse- 
gurada desde que se tire à sorte a(s) fracção(s) 
escolhida (s) para amostra. 

Quanto à fidelidade pode dizer-se que: 


(i) — o erro fundamental depende sobretudo da 
massa da amostra e, em menor medida, 
da massa do lote, 

(ii) — o erro de segregação depende: 


a) da heterogeneidade de distribuição (£). 
Reduz-se homogeneizando o lote, o que 
é importante quando este é uma polpa 
de sólidos finos em suspensão em água; 

bj) da dimensão dos incrementos (+) de- 
positados em cada alvéolo em cada ro- 
tação. Reduz-se aumentando o número 
de alvéolos ou a velocidade de rotação 
ou diminuindo o caudal, 


Conclusão: os divisores rotativos são excelentes, 
recomendando-se o seu uso para amostragem de pol- 
pas e areias secas. Devem evitar-se na amostragem 
de lotes de pequena massa (D, < 1 min), de produtos 
húmidos ou pegajosos, de minérios demasiado gros- 
seiros ou de polpas demasiado instáveis sem prévia 
agitação. 


6. RESOLUÇÃO DOS PROBLEMAS DE AMOSTRAGEM 


6.1 — Introdução 


Em face de cada caso concreto e particular é 
necessário saber como se distribuem os sete erros 


que analisáâmos nos Capítulos precedentes: ED, EL, 
EF, ES, El, EC e EP, quais as precauções a tomar 
para que a amostragem seja justa (pelo menos apro- 
ximadamente), quais os valores a atribuir aos pará- 
metros de funcionamento para que a amostragem 
seja fiel, numa palavra, quais as faltas a não cometer 
para obtenção duma amostra representativa, 

Neste Capítulo, que não é de modo algum um 
receituário, nem tem essa intenção, resolveremos 
sinteticamente um certo número de problemas-tipo 
do lotes manipuláveis e não manipuláveis de várias 
dimensões('). Ver-se-á, inclusivé, que alguns proble- 
mas são insolúveis. 

Diz-se que um problema de amostragem é solú- 
vel quando é possível definir um modo operatório de 
colheita da amostra e atribuir-lhe uma certa repre- 
sentatividade caracterizada por m(ET?) < R2. 


6.2 — Amostragem de lotes sólidos com três dimensões 


Teoricamente será necessário colher incrementos 
distribuídos uniformemente no seio do lote L,. 

Praticamente, dada a natureza não consolidada 
do material a granel(?), não hã qualquer solução 
exceptuando a sua transformação em: 


(i) — lote a duas dimensões com utilização de 
«buldozer»; 


Fig. 13 


(1) Pode afirmar-se que a amostragem dos lotes não manipuláveis é tanto mais fácil e barata quanto menor for o 


número de dimensões do seu modelo. 


(2) As vibrações, choques, pressão e acção perfuradora da sonda tiram qualquer representatividade ao testemunho 


colhido, ao contrário do que sucede com os jazigos mineiros, 
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(iij —lote com zero dimensões (retoma total 
com pá-mecáânica); 


(iii) — lote com uma dimensão temporal (retoma 
e colocação sobre correia transportadora). 


6.3 — Amostragem de lotes sólidos com duas dimensões 


Praticamente é realizável, embora difícil e cara, 
seleccionando Q valores de X e Y segundo uma lei 
de densidade de probabilidade uniforme e colher, em 
cada ponto, um incremento constituído por uma co- 
luna de secção constante, de forma qualquer (por 
exemplo: quadrada) e de altura igual à espessura do 
lote. 


6.4 — Amostragem de lotes sólidos unidimensionais 


É a mais fácil e barata. Consiste na selecção de 
pontos Xp segundo uma lei de probabilidade uniforme, 
ao longo do eixo do lote e colher, em cada ponto, um 
incremento constituído por uma fatia de espessura 
constante e cortando toda a secção do lote, 


6.5 — Concepção e cálculo de uma amostragem de lotes 
com uma dimensão temporal 


6.5.1 — Introdução 


Quando se concebe uma instalação de amostra- 
gem, devem fixar-se certos objectivos dentre os quais 
destacamos os mais importantes: 


1 -— definição do destino e características das 
amostras laboratoriais a obter; 


2 — definição dum nível de representatividade R 2. 
A amostra será representativa se o erro to- 
tal de amostragem ET=(a -a, )/a,, de mo- 
mentos m(ET) e gX(ET), satisfazer a; 


R2=m(ET?) = [m(ET)R+cHET) < R 2 


A representatividade é o único objectivo realista, 
porquanto a precisão, definida por: 


m(ET)=o- (justa) 
c( ET) < R 2 (fiel) 


é impossivel porquanto m(ET) -£o sempre. 
Ora, nós vimos que a correcção da amostragem 


assegura uma justica aproximada, definida por: 


m(ET) « c(ET). 


Vidé: (1) — 3.5. 
(2) — 3.6. 
3) — 3.7 e 4.5, 
(4) — 3,7 e 4,3. 
(5) — 3.8. 
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Como se sabe definir as condições de realização 
duma mostragem correcta, assegurar-se-ã a repre- 
sentatividade através: 


* Correcção: m(ET) = o 
* Fidelidade: R2=9?(ET) < R 2 


6.5.2 — Condições que asseguram a correcção 


1 — Escolha duma lei de integração conveniente, 
assegurando uma densidade de probabilidade 
de colheita uniforme no espaço de L e nula 
no seu exterior (vidé 3.1 e 3.3). 


2 — Escolha duma «découpe» correcta, assegu- 
rando, a todos os fragmentos de L, a mesma 
probabilidade de serem colhidos (vidé 4.5). 


3 — Escolha duma toma correcta (vidé 4.6). 


4 — Anulação dos erros de preparação. 


Nestas condições EC= EP= EZ= O 
o que conduz a: 


m(ET) = O 
(1) (2) (3) 
c(ET)=c(ED) +eoMEI,) +oEF) + 


(4) (5) 
ES) +oUEL,) 


6.5.3 — Condições que asseguram a Fidelidade 


6.5.3.1 — Crédito do erro 


O nível de fidelidade R2 é a variância global 
tolerada na amostragem, no sentido lato. Será pois 
necessário repartir esta variância por cada um dos 
andares de amostragem. Trata-se pois dum problema 
económico cuja solução optimizante será calculada 
para cada caso concreto. 

Como linhas gerais, pode contudo afirmar-se que 
uma amostragem, no sentido restrito, é tanto mais 
cara, quanto maiores forem o caudal, a dimensão dos 
fragmentos e a fidelidade, o que significa que, em 
geral, a amostragem primária será mais cara que 
a secundária, etc. 

Unicamente a título indicativo, baseado na expe- 
riência, poder-se-á dizer que a repartição optimizante 
das variâncias toma em geral os seguintes valores: 


primária: 50 a 75 % de R 2 
secundária: 25 a 15 % de R 2 
terciária: lô5a 5% de R2 


total dos outros andares: 10 a 5% deR2 
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Tendo atribuido a cada andar a variância tole- 
rada, será necessário definir as características e coa- 
dições de funcionamento de cada amostrador, 


6.5.3.2 — Escolha das caracteristicas do amostrador 


Admita-se tratar-se dum amostrador «transver- 
sal» de movimento rectilíineo e seja: 


W — largura da abertura rectangular do 
amostrador [em] 


V — velocidde do amostrador [em .s-1] 

u(t) — caudal no instante t [g.s-1] 

M, —— massa do incremento efectuado no 
instante t [g] 


Me=p(t) Wy/V 
A condição de toma correcta implica 


W>VW, 
V<V, W/W, 


em que: 


Wç=8 d(id=>03em) 
W.=lem (d <0,3cm) 
V. = 40cem.s) 


Quando a toma é correcta, tem-se pois: 
M, > Mag ult) Wo Vo 0.075d ult) 
Verifica-se que não hã interesse em escolher uma 
massa M,. > Mao porqunto: 


* o custo da amostragem a jusante cresce 
com a massa da amostra; 


* para uma mesma massa da amostra M,, 
c2(EI) e g2(ES) são tanto menores quan- 
to menor for a massa do incremento Mo: 


* s2(EF) só depende da massa total da 
amostra. 


Logo, a condição de toma correcta define uma 
massa optimizante do incremento, em cada andar de 
amostragem, que é igual a 


Mgço = 0.075 d p(t)— unidades C, G. S. 


Habitualmente escolher-se-ã W=W, e V=VAD. 


6.5.3.3 — Escolha duma lei de integração 


Se se temem vriações periódicas de qualidade 
escolhe-se uma lei aleatória estratificada. 


Caso contrário escolhe-se uma lei sistemática. 
Em ambos os casos f terá que ser um sub-múl- 
típlo inteiro da duração D do deslocamento do lote. 


6.5.3.4 — Realização da condição de Fidelidade. Valor a 


dar ao intervalo (constante ou médio) É 
Para cada andar n: 


q AET)=o (ED) +o HEI) +o (EL + 
+o EF) +o (ES) 


1. Erro de ponderação ED, 
É imperativo regularizar o caudal e, nes- 
te caso, q (ED) é desprezável. 


2. Erro de integração EL. 

Se se dispõe dos resultados duma expe- 
riência variográfica, calcula-se o *(EI,) para 
diferentes valores de 6 e iterativamente cal- 
cula-se o valor de 8 compatível com o erro 
total, 

Se se não dispõe de informação variográ- 
fica, Gy aconselha, em face da sua experiên- 
cia, que se dê a 6 um valor menor ou igual 
a 10 minutos. 

Para 6 < 10min, o (EL) é desprezável., 


3. Erro de integração do termo periódico EL. 
Se se dispõe dos resultados da variogra- 
fia, determina-se, para a lei escolhida, o va- 
lor máximo tomado por q HEI). 
Se se não dispõe de informação, escolhe- 
-se uma lei aleatória estratificada com uma 
duração de estrato 6 de 5 a l0min. Nestas 
condições, a experiência de Gy mostra que 
q, (El,) é desprezável. 


4, Erro fundamental EF e Erro de Segregacção 
ES. 

Sabemos que q (EF) +o ES)=g HEL). 
Quando se dispõe do variograma, calcula-se 
o, (EL) que é proporcional a 6. 

Quando não se dispõe do variograma, a 
experiência de Gy indica que, em geral, 
(ES) < o EF). 

Donde: q 2(EI,) < 2 c MEF). 

Assim, quando 8 <10min. e o caudal é 
regularizado, a condição de fidelidade é: 


o, (EF) < R,2/2 


Teremos: 
o (EF) = d 3/Me. a R 2/2 


o que conduz a uma massa da amostra no 
andar n 
M,. >2 Cd3R? 


(1) — HA uma excepção, quando se utiliza um amostrador tipo «Godet» e não h& altura disponível para que o «Godet+ 
tenha capacidade necessária para o incremento. Então aumenta-se W (para aumentar a capacidade, para a altura disponi- 


V W 
vel) e Y por forma que — = ——. 
V W 
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a um número de incrementos 
q. 22 C d 3/R/2 Ma, 
em que M; =0.075.d. (E), donde 


q 2 26 C d2/R Ha (E) 


e a um intervalo o. 


9,=D,/Q, < D, KR Htt)/26 C d,2. 


nm 


Se o intervalo assim calculado for inferior a 
10 min., arredondar-se-à para o sub-múltiplo 
de D, mais próximo. 

Se for maior que 10min., tomar-se-á 
como solução o sub-múltiplo de D, mais pró- 
ximo de 10 min. 


5. Exemplo de aplicação, 


Consideremos o minério de ferro do Sul 
de Angola referido em 4.3.6.2. 


Tinhamos estimado d=4,22 em 
C=0.078 g/cm'. 


Suponha-se que se pretende obter uma 
amostra dum período de duas horas (D = 
=7200s) durante o qual o caudal é de 1800 

t/hora (u(t)=0.5X108g.s-1) e que a variân- 
cia relativa tolerada neste andar é R 2=4.10-5. 


Teremos: 


Mo, — 0.075 4.22xX 05x 108 0.16 X 1058 €— 160 kg' 
Q,>26X 0.078X 4.222/4X 10-5x0.5x106=18.07 — 18 


9,=D /Q, =7200/18-400 s=6 min 40 seg. 


Como 8 < 10min é esta a solução. 


Nesta hipótese, a amostra será constituí- 
da por 18 incrementos de 160kg, pesando 
cerca de 2880 kg. 


6.6 — Amostragem de lotes sólidos com zero dimensões 


Consiste na amostragem dum certo número Q de 
unidades inteiras. Desde que se assegure a cada uni- 
dade a mesma probabilidade de ser amostrada (sis- 
temática ou aleatória se não se conhecer N), a amos- 
tragem é correcta. 


On 


Se se conhecer a variância de dispersão dos teo- 
res do lote, a variância da amostra será: 


a = (/Q—1/N) Sur 


Caso contrário não é possível estimar 92,.(1). 


6.7 — Amostragem de lotes manipuláveis 


Já foi tratado no Capítulo 5. 
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O sector siderúrgico nos países do Terceiro Mundo 


(Continuação da pdg. 138) 


Se bem que representando actual- 
mente 70% da população mundial, 
os países do Terceiro Mundo con- 
tribuem apenas com 8% para os 
TOO milhões de toneladas da produ- 
cão de aco bruto, da qual conso- 
mem entre 10 a 13 %, consumo que 
deverá atingir os 29 % até ao ano 
2000. Por outras palavras, os paí- 
ses em questão deveriam criar no- 
vas capacidades da ordem dos 400 
milhões de toneladas, a fim de po- 
derem satisfazer as suas necessi- 
dades sem terem de recorrer à im- 
portação, esta hipótese pode, toda- 
via, considerar-se impraticável na 
medida em que há que distinguir 
entre os países com possibilidades, 
neste domínio, de se lançar neste 
tipo de empreendimentos e aqueles 
para quem tal iniciativa é impra- 
ticável. 

Deste modo, podem eliminar-se 
logo à partida as nações pouco 
povoadas, dado os elevados mon- 
tantes do investimento, os proble- 
mas de matérias-primas, de emer- 
gia e de mão-de-obra que o pro- 
jecto acarretaria, Cite-se, a título 
de exemplo, o Gabão, cujo mer- 
cado ficaria saturado com uma 
produção anual de 200 mil tonela- 
das de aço, sendo 30 mil de re- 
dondos de betão, não se justifican- 
do nestas condições a montagem 
de uma unidade de produção com 
capacidade tão reduzida. Na reali- 
dade, à excepção da Nigéria, não 
parece viável o desenvolvimento 
das actividades siderúrgicas nos 
países da Africa Negra onde, por 
vezes, se revela muito mais ade- 
quada a instalação de pequenas 
unidades de laminagem. 
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e nos países industrializados 


Em contrapartida, a questão 
põe-se de forma diferente em rela- 
ção aos países árabes, os quais, 
embora pobres em minério de ferro, 
dispõem de grandes recursos ener- 
géticos, já se justificando plena- 
mente a instalação de uma indús- 
tria siderúrgica com base no pro- 
cesso de redução directa utilizando 
o gás natural; é o caso dos países 
da OPEP —-Traque, Qatar, Líbia, 
Arábia Saudita e Argélia — que, 
se bem que dispondo de número 
reduzido de consumidores poten- 
ciais, poderão virar-se para a ex- 
portação. 

Não obstante, o mundo árabe, à 
excepção da Argélia onde a pro- 
cura é importante dadas as dimen- 
sões da sua indústria a jusante, 
não entra na categoria dos futuros 
grandes produtores de aço, a qual 
se refere sobretudo aos países com 
mais de cem milhões de habitan- 
tes, que, com um importante mer- 
cado interno a satisfazer, têm 
assegurado o futuro do seu sector 
siderúrgico. É o caso do Brasil, 
india e China que parece disposta 
a tornar-se num grande país side- 
rúrgico (se bem que não caiba na 
classificação nem de país do Ter- 
ceiro Mundo nem de um país em 
vias de desenvolvimento). Efecti- 
vamente, os seus projectos apre- 
sentam tal envergadura que, a con- 
cretizarem-se, modificarão a ima- 
gem da siderurgia mundial. 

A China, quinto produtor mun- 
dial (23,4 milhões de tomeladas) 
entre a RFA e a Itália, é o país 
que reúne maior número de condi- 
ções de desenvolvimento das indús- 


trias siderúrgicas. Com importan- 
tes reservas de minério 1, carvão e 
petróleo, os chineses projectam au- 
mentar a sua produção siderúrgica, 
até ao ano 2000, a um ritmo de 4 
a 6 milhões de toneladas suplemen- 
tares anualmente. Deste modo, e 
através da renovação das unidades 
de produção actuais e construção 
de novas unidades, as capacidades 
de produção passariam de 30 mi- 
lhões de toneladas actualmente 
para 60 milhões em 1985, devendo 
atingir os 88-132 milhões nos pró- 
ximos 22 anos. 

O Brasil ocupa o décimo quarto 
lugar na produção siderúrgica 
mundial (entre a Bélgica e a Es- 
panha) e conta-se igualmente en- 
tre os futuros grandes produtores 
de aço, sendo a exiguidade de pe- 
tróleo compensada pela abundância 
de minério de ferro, que poderá 
utilizar como moeda de troca de 
recursos energéticos. Daí o surto 
da indústria do aco cuja produção, 
2 milhões de toneladas em 1960, 
atingia em 1974 os T milhões, eu 
seja um aumento de 229 %, segui- 
do de um outro de 50 % entre 1974 
e 1977 (11 milhões), num período 
em que os países industrializados 
se viam obrigados a reduzir as 
suas capacidades. Entretanto, pre- 
vê-se ainda que sejam realizados, 
até 1981, investimentos no montan- 
te de 3,5 mil milhões de dólares 
tendo em vista duplicar a produção 
num periodo de 10 amos. 


1 A sua qualidade não é das melho- 
res; entretanto, serão exploradas no- 
vas minas. As principais fontes de 
abastecimento da China no exterior 
são a Austrália e o Brasil. 
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A índia, décimo sexto produtor 
mundial (10 milhões de toneladas 
— entre a Espanha e a Austria), 
pertence também à categoria dos 
futuros grandes produtores de aço; 
todavia, terá de percorrer ainda 
um longo caminho uma vez que 
o consumo local atinge apenas 
16 kg/ano/habitante, contra 100 kg 
no Brasil, 370 em França, 512 no 
Japão e 618 nos Estados Unidos. 
O país pretende aumentar a sua 
capacidade total para 60-75 mi- 
lhões de toneladas anuais no ano 
2000, 

Face ao exposto, e independen- 
temente da importância dos pro- 
jectos previstos para a América 
Latina, Extremo Oriente” e outros 
países, conclui-se que está longe de 
ser possível atingir os 400 mi- 
lhões de toneladas sunlementares 
que permitiriam aos países do Ter- 
ceiro Mundo atingirem o nível de 
auto-abastecimento. Aliás, do pon- 
to de vista dos peritos, e excepto 
no que diz respeito à China, a reali- 
dade apontaria para valores me- 
nores em cerca de 50% aos esti- 
mados nas previsões referidas, 
Com efeito, só o= países industria- 
lizados têm possibilidades, do pon- 
to de vista técnico, de construir 
rapidamente tão importantes capa- 
cidades, já que os países em vias 
de desenvolvimento não têm condi- 
cões, inclusivamente mo aspecto fi- 
nanceiro, para, por si sós, levarem 
a cabo tal empreendimento. 

Cada país tem os seus proble- 
mas próprios em função dos res- 
pectivos mercados, condições dos 
povos e recursos em matérias-pri- 
mas, verificando-se que, na maior 
parte dos casos, se tem posto em 
funcionamento unidades de produ- 
ção sem quadros técnicos e operá- 
rios especializados, o que dá ori- 
gem a uma série de problemas 
dificeis de solucionar, E o caso do 
Qatar (200 mil habitantes) e da 
minifábrica de UMM Said, que en- 
trou em servico este ano e cujo 
funcionamento será assegurado, 
por um periodo de oito anos, pelo 
grupo japonês Kobe Steel, respon- 
sável pela sua construção. Durante 
os trabalhos, a Kobe forneceu uma 
equipa composta por mais de 100 
pessoas encarregada de vigiar vá- 
rias centenas de operários vindos 
da Coreia do Sul, Filipinas, Índia 
e de outros países. Os engenheiros, 
contramestres e operários qualifi- 
cados foram recrutados no Egipto 
e Extremo Oriente. 

A nova geografia do aço encon- 
tra-se, pois, longe de constituir 
uma ameaça para os países indus- 
trializados. Nos países de sector 
siderúrgico recente, as necessida- 
des tendem a desenvolver-se a 
riimo mais acelerado que o do 
aumento da capacidade de produ- 
ção, cujo arranque é, por outro 
lado, muito mais moroso; nenhum 
deles estará, pois, à altura de fazer 
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concorrência aos Estados Unidos, 
Japão ou à CEE em matéria de 
tonelagem. Estes países reduzirão, 
entretanto, as suas importações, o 
que irá afectar consequentemente 
as empresas do ramo japonesas e 
do Mercado Comum, exportadores 
tradicionais. Para estes dois blo- 
cos subsiste ainda um factor mais 
inquietante, o de verem um dia os 
países do Terceiro Mundo cortar 
o fornecimento de minério com o 
objectivo de fazer valer as vanta- 
gens de que dispõem. Refira-se, a 
título de exemplo, que, no caso de 
o Brasil decidir bruscamente in- 
terromper as suas exportações (60 
milhões de toneladas/ano, num to- 
tal de 300 milhões comercializados 
em todo o mundo), a Europa — e 
o Japão —, que importa 130 mi- 
lhões de toneladas, ficaria a bracos 
com a necessidade imperiosa de 
encontrar nova fonte de forneci- 
mento que para si camalizazse 05 
30 milhões de toneladas de que 
não poderia prescindir. 


Produção mundial de aço por regiões 
(Valores de 1977 em toneladas) 


Europa Ocidental 155 210 000 
Europa de Leste 147 000 000 
América do Norte 126 810 000 
Oceânia 7500 000 
América Latina 21 870 000 
dos quais América 
Central 15 000 
Argentina 2 700 000 
Brasil 11 200 000 
Chile 560 000 
Colômbia 400 000 
Cuba SO OO 
México 5 900 000 
Peru 320 000 
Uruguai 20 000 
Venezuela 814 000 
África 8 400 000 
dos quais Africa 
do Sul 7300 000 
Argélia 400 000 
Rodésia 500 000 
Tunísia 100 000 
Outros 150 000 
Ásia 146 600 000 
dos quais China 23 400 000 
Coreia do Norte 3 200 000 
Coreia do Sul 4 200 000 
Filipinas S00 000 
Formosa 7100000 
Hong-Kong' 120 000 
índia 10 000 000 
Indonésia 100 000 
Japão 102 400 000 
Malásia 250 000 
Singapura 200 000 
Tailândia 300 000 
Outros 350 000 


Isto significa que seriam neces- 
sários vários anos de trabalho de 
abertura de minas para o conse- 
guir. Daí os interesses no domínio 
das indústrias extractivas adquiri- 
das, no decurso dos anos, pelas 


siderurgias da CEE nos países pro- 
dutores, onde controlam actual- 
mente uma produção da ordem dos 
17 milhões de toneladas. 

Nem só a siderurgia europeia 
foi afectada pela crise que tem 
grassado nos últimos quatro anos 
neste sector; o Japão, que contri- 
buiu para o agravamento da crise 
na Europa, sentiu também, ainda 
que menos acentuadamente, os 
seus efeitos. Assim, a produção si- 
derúrgica japonesa foi reduzida de 
apenas 10% no período 1973-T6, 
tendo o consumo local de aco bai- 
xado de 27,5 %. 

A produtividade na indústria si- 
derúrgica nipónica mantém-se im- 
batível, embora inferior à até aqui 
verificada, tendo atingido em 1977 
as 327 toneladas de aço anuais por 
trabalhador, contra 361 em 1973. 
Em contrapartida, a produtividade 
da CEE não ultrapassou, em mé- 
dia, 175 toneladas em 1977, situan- 
do-se a Franca em último lugar 
com 149 toneladas. 

A fim de sobreviver à queda do 
consumo e à guerra dos preços, as 
empresas do sector recorreram a 
Bruxelas — cuja resposta consis- 
tiu no Plano Davignon que, aliás, 
todas elas procuram ludibriar —, 
endividaram-se e recorreram ao 
despedimento de trahalhadores. O 
endividamento do sector atingiu 
em França, no ano de 1977, um 
recorde de 111% do seu volume 
de negócios (contra 87 % um ano 
antes, data em que correspondia, 
em Ttália, uma taxa de 104 & con- 
tra apenas 22 e 21%, respectiva- 
mente na República Federal da 
Alemanha e nos Estados Unidos). 
Quanto à redução de efectivos, só 
em 1977 a França começou a re- 
correr a este tipo de soluções (ten- 
do procedido a 11500 despedimen- 
tos), uma vez que até então o Es- 
tado tudo fez para evitar tal me- 
dida. Entretanto, para fazer face 
à crise, a maior parte dos restan- 
tes países europeus não hesitou 
em fechar os olhos a problemas 
sociais. No periodo 1974-77T, as si- 
derurgias alemãs despediram 25 
mil efectivos muito embora tives- 
sem procedido à reestruturação do 
seu aparelho de produção antes da 
crise, que se traduzira por um9 
contracção de efectivos da ordem 
dos 20 mil entre 1960 e 1974 (con- 
tra 500 em França). Esta situação 
leva o banco privado alemão 'Trim- 
kaus und Burkhardt a concluir, 
num estudo recentemente elabora- 
do sobre a situação do sector side- 
rúrgico no país, que o período mais 
crítico teria terminado, 


2 As capacidades atingiriam, até 1985, 
4,5 milhões de toneladas na Argentina, 
16 milhões no México, 10 milhões na 
Venezuela e 20 milhões na Coreia do 
Sul (perspectivas que parecem merece- 
doras das maiores reservas). 


fContinua na pdg." 208) 
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Representatividade da amostragem de jazidas minerais 
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RESUMO 


São analisados aspectos da representitividade da amostragem de jazidas minerais destinada a ensaios 
mineralúrgicos e metalúrgicos. 

E definida população de referência e apresentada metodologia para a reproduzir por amostragem 
e apresenta-se um exemplo real extraído do estudo de extensa bacia carbonifera Brasileira. 


ABSTRACT 


Comments are developed on the representativity of ore bodies sampling for obtaining data concerning 
the behaviour on concentration and metallurgical processes. 

The reference population is defined and some methods for reproducing it are presented. 

An example taken from a study carried out in « Brazilian coal basin is summarized, 


1. INTRODUÇÃO 


O problema que aqui se aborda não diz respeito ao reconhecimento de uma jazida por meio da amos- 
tragem, estudo este efectuado nas fases de pesquisa e reconhecimento. 

Supõe-se, pelo contrário, que esse reconhecimento está feito (com maior ou menor rigor), que se defi- 
niram, com certos erros, pelo menos, as dimensões e forna da jazida, o teor médio global, a massa mine- 
ralizada global, isto é, se avaliaram as reservas; pode ainda haver-se mesmo estabelecido determinada 
tipologia, baseada, por exemplo, nos teores tal-qual e/ou na mineralogia e/ou mineralurgia, ou ainda em 
propriedades mecânicas, etc. 

De posse de todos estes dados, trata-se sim do problema de colher amostras da jazida ou de parte 
importante dela, destinadas a estudos mineralúrgicos que, como se sabe, se iniciam por testes laboratoriais 
descontínuos (testes de bancada) para se passar, seguidamente, aos testes piloto ou semi-industriais contínuos, 
a estes dois escalões de ensaios correspondendo, como é óbvio, volumes de amostra muito diferentes. 


2. POPULAÇÃO DE REFERÊNCIA 


Na realidade, a pesquisa e reconhecimento forneceram uma população de referência que resultou da 
amostragem já efectuada na jazida. Mediante certo modelo de função aleatória compativel com os dados e útil 
para orientar a inferência estatística terá sido possível, na maioria dos casos, construir modelo(s) de vario- 
grama(s) baseado(s) nos variogramas experimentais e, mesmo (se necessário) estimar uma deriva e reali- 
zar, em qualquer caso, kKrigagem da jazida, tomando por base descritização relacionada com a unidade ou 
o bloco tecnológico, isto é, aquele volume sobre o qual se poderão tomar opções nos actos da futura explo- 
ração, 

Consequentemente, referido(s) a tal bloco tecnológico, dispõe-se de histograma(s) que são, evidente- 
mente, a tradução da população dos valores krigados, a que corresponde certo valor médio m=E(Z ) = 
= E(Z) e certa variância de dispersão D(Z0). 


É claro que se tiver havido necessidade ou conveniência em definir, já nesta fase, uma tipologia, a 
jazida terá sido dividida em sub-jazidas, a cada uma das quais corresponderá um histograma para o atri- 
buto considerado relevante, 
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Além disso. poder-se-á estar interessado em mais de um atributo. Por exemplo, no caso dos carvões, 
poderá interessar os teores de cinzas, matérias voláteis, carbono fixo, enxoire, etc. 

E todos estes diferentes atributos poderão ser necessários para caracterizar qualquer dos minérios 
que a tipologia definiu na jazida. Dispor-se-á, assim, às vezes, de colecção relativamente numerosa de histo- 
gramas: uns tantos para cada tipo de minério. 


3. OBJECTIVOS DA AMOSTRAGEM. SUA CONSTITUIÇÃO 
3.1. Primeiro objectivo: Investigação do tratamento mais adequado 


A representatividade de uma amostragem não tem significado positivo se não se houver previamente 
definido o seu objectivo. 
Suponha-se então que tal objectivo se pode genericamente caracterizar do seguinte modo: 


«Pretende-se obter amostragem para investigar, em primeira fase, da adequação de determinado(s) 
processo(s) mineralúrgico(s) e/ou metalúrgicos(s)». 

Trata-se, obviamente, de realizar ensaios de bancada diversos, os quais deverão estudar a resposta, 
ao tratamento em questão, do espectro que o histograma reflecte. 

Quer dizer, o histograma traduz uma sub-tipologia, expressa pelo parâmetro respectivo, a qual terá 
de ser ensaiada em toda a extensão. 

Há então que investigar dominio(s) A (Zan Zmax) Por onde se espraia(m) o(s) atributo(s) julgado(s) 
relevante(s). Supondo, por simplicidade, que um único atributo é relevante, deverá, então, dividir-se o domi- 
nio respectivo, em intervalos disjuntos AZ,, tais que |) Pa; = A, e escolher, na jazida, locais onde haja for- 
tes probabilidades das amostras aí colhidas exibirem valores do atributo que caiam dentro de cada um dos 
intervalos AD, assim arbitrariamente escolhidos. 

A tais amostras chamaremos incrementos. Os respectivos pesos variarão, conforme os casos, entre 05 
e 150 Kg. 


Dois casos se poderão então dar: 


1.º — Há boa acessibilidade à jazida: por exemplo, os testemunhos das sondagens existem, têm largo diá- 
metro, estão ordenados e bem referenciados, não se encontram alterados e será, portanto, possível, 
fazer colheita de amostras nos próprios testemunhos (esta observação chama a atenção para a necessi- 
dade de projectar convenientemente as campanhas de sondagens, posteriores às da pesquisa inicial, 
tomando em consideração os diâmetros necessários para virem a fornecer os incrementos referidos). 
Ou ainda, outro exemplo, a jazida dispõe de ossatura básica de preparação para desmonte — galerias, 
travessas, chaminés, devidamente amostradas. Será, então, possível, colher aí tantos incrementos quantos 
os necessários para satisfazer o desideratum — uma amostra. de teor AZ, representando cada inter- 
valo A Z. | 

Será desejável que a amplitude dos intervalos AZ, não seja pequena, por razões óbvias de eco- 
nomia; por outro lado, deverá constituir-se a amostra representativa de cada intervalo AZ, com incre- 
mentos colhidos em diferentes locais, previamente seleccionados nas carotes ou nas galerias, travessas, 
chaminés, etc., para que atributos ocultos, mas existentes, venham a ser considerados na investigação, 
Teoricamente, cada amostra representativa de AZ, deverá ser constituída por k incrementos de teor 
médio m, que satisfaçam, por exemplo, às condições 


a! ( My — A Z, ) —+» () Dad = 
al ( Mj— & Z) —» () a! > 0 
Mas aqui surgem certas dificuldades, a saber: 


(ij) Cada incremento será constituido por conjunto de peças minerais a granel, de diferentes cali- 
bres. Deverão, assim, por si próprios, satisfazer às condições estabelecidas pela teoria da amos- 
tragem de Gy, isto é, deverão ter pesos tais que os tornem representativos das granulometrias 
que eles próprios exibem. 


Antes de mais, deve considerar-se duas situações favoráveis: 


— O processo de tratamento que se pretende ensaiar exige, a priori, que se trabalhe com gra- 
nulometria fina (por exemplo, inferior a 1 mm, ou 1004 ou 504, etc.) e há disponíveis teste- 
munhos de sondagens. Neste caso, a simples extracção de um segmento de testemunho e a sua 
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